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FOREWORD 


The  1998  version  of  the  aeroprediction  code  (AP98)  was  based  on  the  use  of  slender 
body  and  linearized  theories  to  predict  aerodynamics  at  low  angles  of  attack.  A  missile 
component  wind  tunnel  data  base  taken  in  the  1970s  at  NASA  Langley  Research  Center  was 
used  to  predict  aerodynamic  nonlinearities  as  a  function  of  Mach  number,  angle  of  attack,  and 
various  missile  geometric  parameters.  A  more  recent  data  base  also  conducted  at  NASA  was 
made  available.  This  data  base  was  therefore  used  to  refine  the  empirical  constants  used  in  the 
aeroprediction  code  to  predict  the  aerodynamic  nonlinearities.  This  report  documents  the  new 
empirical  constants  derived  and  the  improvements  in  accuracy  of  normal  force  coefficient 
afforded  by  this  new  set  of  constants. 

The  work  described  in  this  report  was  supported  through  the  Office  of  Naval  Research 
through  the  Surface  Weapons  Systems  Technology  Program  managed  at  the  Naval  Surface 
Warfare  Center,  Dahlgren  Division  (NSWCDD)  by  Mr.  Robin  Staton.  Tasking  from  this 
program  was  provided  by  Mr.  Roger  Horman  and  Mr.  John  Fraysse.  Also,  some  support  was 
provided  by  the  Marine  Coips  Weaponry  Technology  Program  managed  at  NSWCDD  by 
Mr.  Craig  Melton.  The  authors  express  appreciation  for  support  received  in  this  work. 


Appreciation  is  also  expressed  to  Mr.  Tom  Hymer  (G23)  who  assisted  the  lead  author  in 
the  last  week  of  the  writing  of  this  report  to  fine-tune  some  of  the  empirical  constants. 
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1.0  INTRODUCTION 


The  1998  version  of  the  NSWCDD  Aeroprediction  Code  (AP98)1  is  the  most  complete 
and  comprehensive  semiempirical  code  produced  to  date.  It  includes  the  capability  to  predict 
planar  aerodynamics  in  the  roll  positions  of  =  0  deg  (fins  in  ‘  +  or  plus  orientation  as  viewed 
from  the  rear  of  the  missile)  and  0  =  45  deg  (fins  in  “x”  or  cross  roll  orientation  as  viewed  from 
the  rear  of  the  missile)  over  a  broad  range  of  flight  conditions  and  configuration  geometries  with 
good  average  accuracy,  computational  times  and  ease  of  use.  Flight  conditions  include  angles  of 
attack  (AO A)  up  to  90  deg,  control  deflections  of  up  to  ±30  deg,  and  Mach  numbers  up  to  20. 
Configuration  geometries  include  axisymmetric  and  nonaxisymmetric  body  shapes  with  sharp, 
blunt,  or  truncated  nose  tips,  with  or  without  a  boattail  or  flare.  Up  to  two  sets  of  planar  or 
cruciform  fins  are  allowed.  New  technology  has  recently  been  developed2  to  allow  both  six-  and 
eight-fin  options  in  the  fin  considerations  as  well.  Average  accuracies  are  ±10  percent  for 
normal  and  axial  force  and  ±4  percent  of  body  length  for  center  of  pressure.  By  average 
accuracy  is  meant  that  enough  AOAs  or  Mach  numbers  are  considered  to  get  a  good  statistical 
sample.  On  occasion  a  single  data  point  can  exceed  these  average  accuracy  values.  Ease  of  use 
has  been  significantly  enhanced  over  older  versions  of  the  Aeroprediction  Code  (APC)  through  a 
personal-computer-based  pre-  and  post-processor  package.3  This  package  has  allowed  inputs  for 
configuration  geometries  to  be  simplified  significantly  by  many  automated  nose  shape  options. 


While  the  AP98  is  a  very  powerful  tool,  several  limitations  and  areas  of  improvement 
still  remain.  Most  of  these  needs  are  driven  by  the  desire  of  future  weapon  designers  to  perform 
trade  studies  on  new  and  innovative  concepts  that  may  fall  outside  of  the  current  capability  of  the 
AP98.  An  example  of  this  type  of  requirement  is  the  multi-fin  requirement  that  has  just  been 
completed.2  Another  example  of  this  type  of  requirement  is  to  include  the  capability  to  deflect 
the  rear  segment  of  a  fin  (sometimes  referred  to  as  flaperon  or  aileron)  for  control,  as  opposed  to 
the  entire  fin. 

In  addition  to  configuration  design  flexibility,  there  are  several  areas  where  improvement 
in  the  aerodynamics  computational  process  is  needed.  The  semiempirical  model  for  the  wing 
and  wing-body  interference  aerodynamics  was  based  primarily  on  missile  component  data 
bases4'6  where  the  parameter  r/s  was  a  constant  value  of  0.5  (see  Figure  1  for  nomenclature). 
More  recently,  a  new  missile  component  data  base  has  been  made  available  where  data  was 
measured  for  wing-alone  and  wing-body  configurations  with  r/s  =  0.25,  0.33,  and  0.5.  This  new 
data  base  should  therefore  allow  refinements  in  the  AP98  methodology  for  the  wing  alone  as 
well  as  for  the  effects  due  to  r/s.  Other  refinements  not  explicitly  considered  in  the  AP98  include 
an  approximate  way  to  account  for  internal  shock  interactions  and  improvements  to  the  wing-tail 
interference  model.  The  wing-tail  interference  model  was  also  based  on  a  limited  data  set. 
Issues  such  as  internal  shock  interactions  and  wing-tail  interference  should  be  amenable  to 
treatment  by  computational  fluid  dynamics  (CFD)  codes.9' 


1 


FIGURE  1.  TYPICAL  AXISYMMETRIC  WEAPON  CONFIGURATION  GEOMETRY 
OPTIONS  AND  NOMENCLATURE 


The  purpose  of  this  report  is  to  examine  the  recent  data  base  in  Reference  7  and  make 
appropriate  improvements  in  the  AP98  methodology.  Once  these  improvements  have  been 
made  they  will  be  validated  against  a  broad  range  of  configurations  other  than  those  upon  which 
the  data  was  measured.  Assuming  the  improvements  improve  the  accuracy  of  the  predictions  for 
the  validation  cases,  the  methodology  will  be  integrated  into  the  AP98  and  will  be  transitioned  as 
the  AP02  in  fiscal  year  2002. 


2.0  MODIFICATIONS  TO  THE  AP98  BASED  ON  RECENT  WIND  TUNNEL  DATA 


The  AP98  uses  the  so-called  direct  approach  to  estimating  the  nonlinear  terms  of  each  of 
the  normal  force  and  center  of  pressure  components.  By  direct  approach  is  meant  that  each  of 
the  terms  in  Equation  (1)  below  is  broken  down  into  a  linear  and  a  nonlinear  component.  In 

Equation  (l),12 


+  CN  +  CN 


B(W) 


+  CNt,bi  +CNB(T1  +CNt 


(1) 


it  is  understood  that  the  CN#(V)  term  (not  shown)  is  included  in  the  CNb<ti  term.  As  an  example 
of  the  direct  approach,  consider  the  CNw  bi  component  of  Equation  (1),  which  can  be  expanded 

to 
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CNW(B)  =  [(CN0  )L  +  (CNa  )NL  l  {[(KW(B)  )sbt  +  (^WCB)  )NL  ] 

(2) 

+  (C1  [k  W(B)  IbT  +  C2  ^  T 

VAref 

The  linear  or  small  angle  of  attack  terms  of  Equation  (2)  are  estimated  by  linear  theory  (LT)  or 
slender  body  theory  (SBT).  This  gives  the  APC  a  good  fundamental  basis  for  its  aerodynamic 
estimates.  The  nonlinear  corrections  due  to  higher  AOA  or  control  deflection  are  each  estimated 
directly  from  component  wind  tunnel  data  bases.  Each  of  the  other  terms  in  Equation  (1)  is 
treated  in  a  similar  fashion  to  Equation  (2)  in  the  actual  implementation  into  the  APC. 

The  primary  data  bases  used  to  define  the  nonlinear  terms  of  Equations  (1)  and  (2)  were 
those  of  References  4-6  (see  Figure  2).  In  those  references,  measurements  of  static 
aerodynamics  for  wing  alone,  body  alone,  and  wing  in  close  proximity  to  the  body  were  made  as 
a  function  of  roll  angle,  AOA,  Mach  number,  and  control  deflection.  When  the  wing  in 
proximity  to  the  body  measurements  were  made,  two  sting  measurements  were  simultaneously 
made  to  record  the  load  on  the  wing  close  to  the  body  and  the  load  on  the  body  close  to  the  wing. 
This  allows  direct  measurement  of  CN  ,  CN  ,  and  CN  of  Equation  (1).  Knowing  CNb  and 

CNw  from  previous  wind  tunnel  measurements,  values  of  the  interference  terms 


can  be  obtained.  ACNb(w)  of  Equation  (3B)  is  obtained  from 

ACNB(W)  =CN  -CNb  “CNw,,,, 


(3C) 


Reference  13  performed  a  qualitative  error  analysis  of  obtaining  Kw(b>  and  KB(w>  from 
Equations  (3)  and  References  4—6.  The  conclusions  of  that  study  were: 

a)  KW(B)  estimates  are  reasonable  for  AR  <  1.0  at  all  M  and  a.  More  scatter  occurs  for 
AR  >  1 .0. 

b)  KB(w>  estimates  are  usable  for  lower  AR  and  M.  Best-guess  judgement  is  needed  for 
other  conditions. 


3 


NSWCDD/TR-99/1 16 


Taper 

Ratio 


0.0 


0.25 


Aspect  Ratio 
0.5  1.0  2.0  4.0 


0.5 


z  :  i 


m  m 


1.0 


CD 


18 


4.5  2.25  1.125 


Surface  Area  (in. ) 


- V 

Fixed  Fins 


Control  Fins 
I  !  6f=  -40°  +40° 

+4=1 

Hinge  Lines 


Body 

3  cal.  Tangent 
Ogive  Nose 
12.33  cal.  Length 
3  in.  Diameter 


FIGURE  2A.  MODELS  USED  IN  LANGLEY4  WING-BODY  TESTS 
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FIGURE  2B.  MODELS  USED  IN  STALLINGS  AND  LAMB5  WING-ALONE  TESTS 
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These  conclusions  were  primarily  driven  by  the  small  fins  of  the  Reference  4  data  base,  which 
were  required  in  order  to  keep  r/s  of  Figure  1  constant  at  a  value  of  0.5.  As  a  result. 
Reference  13  recommended  that  any  future  component  test  data  for  engineering  codes  should 
have  reasonably  sized  fins  in  relation  to  the  body  planform  area  and  mounted  in  the  mid-body 
region  to  capture  most  of  the  afterbody  carryover  aerodynamics. 

The  recent  wind  tunnel  test,7  carried  out  jointly  by  then  McDonnell  Douglas  Aerospace 
and  NASA/LRC,  followed  the  recommendation  of  Reference  13  for  larger  wings,  as  most  of  the 
wing  planforms  were  equal  to  or  larger  than  those  of  Reference  4.  Figure  3  shows  the 
configuration  geometry  and  dimensions  of  the  body  and  wing  planforms  tested,  and  Figure  4 
shows  a  scaled  pictorial  view  of  the  Reference  7  and  Reference  4  configurations  tested.  These 
can  be  compared  directly  to  the  test  of  Reference  4  and  Figure  2.  Several  points  are  worthy  of 
note.  First,  the  wind  tunnels  used  in  Reference  4  were  the  supersonic  tunnel  at  NASA/LRC  and 
the  subsonic  facility  at  NASA/AMES,  whereas  the  tunnels  used  in  Reference  7  were  both  at 
NASA/LRC.  The  supersonic  tunnel  was  the  same  as  that  of  Reference  4,  but  the  transonic  and 
subsonic  data  were  from  the  NASA/LRC,  8-ft  transonic  facility.  Second,  the  Mach  number 
range  of  Reference  4  was  from  0.6  to  4.6,  whereas  that  of  Reference  7  was  from  0.6  to  3.95. 
Third,  the  AOA  range  of  Reference  4  was  from  0  to  40,  whereas  the  range  of  Reference  7  was 
from  0  to  25  to  30  deg,  depending  on  the  Mach  number.  The  nose  shape  of  both  bodies  were 
identical,  but  the  Reference  7  tests  had  an  afterbody  length  1  in.  shorter  than  that  of  Reference  4. 
Hence,  Reference  4  data  is  for  a  12.33-caliber  body  and  Reference  7  data  is  for  a  12.0-caliber 
body. 5  The  single  wing  planform  areas  of  the  Reference  4  data  base  that  were  tested  in 
conjunction  with  the  body  varied  in  area  from  1.125  sq.  in.  to  18  sq.  in.,  whereas  those  of 
Reference  7  varied  in  area  from  2.25  sq.  in.  to  20.25  sq.  in.  However,  the  largest  fin  of 
Reference  4  was  for  aspect  ratio  of  0.25  and  the  data  base  was  not  complete.  Hence,  effectively 
the  wing  area  of  the  Reference  4  data  varied  from  1.125  to  9  sq.  in.;  so,  in  effect,  the  wing  sizes 
of  the  Reference  7  data  were  about  twice  the  size  of  those  tested  in  Reference  4.  Another  major 
difference  was  the  fact  that  the  wing-alone  data  used  in  the  AP98  was  based  primarily  on 
Reference  5  and  Reference  14.  The  Reference  5  data  was  taken  using  a  sting  mount  in  the  tunnel 
and  integrated  pressure  data.  While  it  is  believed  the  wing  alone  data  base  of  Reference  5  could 
be  slightly  low  in  some  cases  because  of  thickness  effects,  the  author  believes  this  is  the  best 
wing-alone  data  base  available.  The  Reference  7  data  for  wing-alone  aerodynamics  was  taken 
based  on  the  same  wings  of  Figure  3  that  were  tested  on  the  body  alone,  in  contrast  to  the 
References  4  and  5  data,  where  different  size  wings  were  tested  for  wing  alone  and  wing  in 
conjunction  with  the  body  data  (because  of  requirements  for  many  pressure  taps  in  the  wing- 
alone  measurements).  The  wing-alone  data  of  Reference  7  was  obtained  on  a  splitter  plate, 
versus  a  sting  in  the  Reference  5  data,  and  as  will  be  shown  later,  this  mount  arrangement 
apparently  caused  measurement  errors  at  some  conditions.  Finally,  only  two  fins  were  mounted 
on  the  body  in  the  $  =  0  roll  orientation  in  the  Reference  7  tests,  whereas  four  fins  were  mounted 
on  the  body  in  the  Reference  4  data  base  and  roll  angle  was  varied  as  well. 

Both  data  sets  (References  4  and  5  and  Reference  7)  measured  aerodynamics  for  body 
alone,  wing  alone,  and  wing  in  conjunction  with  the  body.  As  a  result,  the  approach  taken  here 
will  be  to  apply  the  AP98  methodology  directly  to  the  Reference  7  normal  force  coefficients  for 
body  alone,  wing  alone,  and  wing  in  conjunction  with  the  body.  Based  on  comparison  of  the 
AP98  to  the  data,  areas  where  improvements  are  needed  will  be  identified.  Modifications  to  the 
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FIGURE  3A.  FORCE  AND  MOMENT  SIGN  CONVENTIONS  OF  NASA/MDAC  DATA  BASE7 


FIGURE  3B.  WING  PLANFORMS  TESTED  IN  NASA/MDAC  DATA  BASE7 
(ALL  DIMENSIONS  IN  INCHES) 
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FIGURE  4.  SCALED  GEOMETRY  COMPARISON  WITH  TRI-SERVICE  MODEL4 

AND  REFERENCE  7  MODEL 


AP98  methodology  will  then  be  made.  These  modifications  will  then  be  checked  out  on  the 
Reference  4  and  7  data  bases  and  then  on  several  missile  configurations  outside  the  Reference  4 
and  7  data  bases.  Any  fine-tuning  necessary  in  the  methodology  improvements  will  be 
conducted  as  part  of  the  methodology  development  process.  If  the  improvements  prove  to  be 
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effective  in  an  overall  sense,  the  new  methodology  will  be  made  a  part  of  the  next  version  of  the 
APC,  which  is  scheduled  to  be  released  as  the  AP02. 


2. 1  BODY- ALONE  MODIFICATIONS 

The  present  body-alone  static  aerodynamics  are  computed  using  linearized  theories  at 
low  AOA  and  a  modified  version  of  the  Allen  Perkins  viscous  crossflow  for  the  nonlinear  AOA 
aerodynamic  terms.  One  of  the  keys  in  obtaining  accurate  aerodynamics  is  in  obtaining  accurate 
values  of  the  critical  crossflow  Reynolds  number  and  Mach  number.  These  parameters  are  of 
primary  importance  at  low  Mach  number.  For  Mach  number  2.0  and  greater,  they  have  little 
influence  on  the  aerodynamic  terms.  The  AP98  currently  uses  a  value  of  RNc  =  180,000  and 

M  =  0.1  as  standard  values.  However,  the  user  is  allowed  to  change  RN(.andMNc  to 

NC 

specified  values. 

In  comparing  the  AP98  to  the  body-alone  wind  tunnel  data  of  Reference  7,  good 
agreement  in  center  of  pressure  and  normal  force  were  obtained.  Average  errors  of  normal  force 
were  less  than  6  percent  and  center  of  pressure  less  than  !4  caliber  or  2  percent  of  the  body 
length.  These  average  errors  were  calculated  using  optimum  values  of  the  critical  crossflow 
Mach  number  and  Reynolds  number,  which  is  quite  important  for  M«  <1.2  comparisons.  RNc 

was  a  constant  330,000,  and  MNc  varied  from  0  at  M»  =  0.6  to  0.06  at  M»  =  0.9.  Also,  error 

values  were  calculated  at  each  5-deg  AOA  at  all  Mach  numbers  where  data  was  available.  This 
gave  a  total  of  40  data  points,  sufficiently  large  to  get  a  good  statistical  average  error. 

In  viewing  the  individual  comparisons,  it  was  clear  that  a  couple  of  minor  problems 
existed,  which  if  corrected,  could  improve  these  average  errors  somewhat.  The  first  one  has  to 
do  with  the  fact  that  the  current  body-alone  methodology  for  implementing  compressibility 
effects  into  the  nonlinear  normal  force  term  could  be  improved  upon.  The  present  methodology 
for  the  body-alone  aerodynamics  in  the  normal  plane  is: 

CNB=CNL+TiCdcsin2a-^ 

Aref 

(Xcp  )l  CNl  (XCp  )NL  CNnl 


CMb=-CNb(Xcp-X0)  (4C) 

In  addition,  an  empirical  table  of  center  of  pressure  shifts  was  used  for  the  body-alone  to 
partially  account  for  physics  not  adequately  accounted  for  in  the  determination  of  center  of 
pressure.  These  physics  include  the  following:  transonic  flow  where  shock  waves  can  stand  on 
the  body,  the  fact  that  the  linear  theory  center  of  pressure  does  not  stay  constant  as  is  presently 
assumed,'  and  the  fact  that  the  center  of  pressure  moves  in  a  parabolic  fashion  (versus  a  weighted 


(4A) 

(4B) 
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average  as  represented  by  Equation  (4))  from  its  value  at  a  =  0  to  the  centroid  of  the  planform 
area  at  a  high  AOA,  say  45  deg. 


Three  slight  changes  in  the  Reference  1  methodology  are  being  implemented  as  a  result 
of  comparisons  to  the  Reference  7  data  base.  The  first  has  to  do  with  the  value  of  tj.  r|  is  the 
normal  force  of  a  circular  cylinder  of  given  length-to-diameter  ratio  to  that  of  a  circular  cylinder 
of  given  length.  T|o  is  the  value  of  rj  at  Mn  =  0.  At  present, 


n= 

=i 


i-n0 

1.8 


Mn  +T|0forMN  <1.8 


forMN  >1.8 


(5) 


Also,  T)  is  automatically  set  to  one  if  M«  >  2.75.  This  last  condition,  where  r\  is  automatically  set 
to  one,  appears  not  to  be  necessary.  In  other  words.  Equation  (5)  is  allowed  to  be  the  sole 
determination  of  the  value  of  t|.  This  change  mainly  affects  normal  force  results  for  conditions 
just  above  the  cutoff  Mach  number  of  2.75.  Figures  5  and  6  compare  the  revised  methodology  of 
removing  the  r)  =  1  for  M»  >  2.75  condition  and  using  only  Equation  (5)  for  normal  force,  versus 
the  current  AP98  approach.  Figure  5  is  the  Reference  7  data  base  and  Figure  6  is  the  Reference  4 
data  base.  Results  are  shown  only  for  Mach  numbers  above  2.75.  Both  normal  force  coefficient 
and  center  of  pressure  are  given.  Note  that  the  revised  method,  which  will  be  incorporated  as  a 
part  of  the  AP02,  shows  improvement  in  comparison  to  both  the  Reference  7  and  Reference  4 
data  base.  The  average  6  percent  error  of  the  AP98  compared  to  the  Reference  7  data  base  is 
reduced  to  an  average  error  of  about  4  percent.  Also,  some  improvement  in  the  average  error 
comparisons  of  the  Reference  4  data  base  is  obtained,  although  this  error  was  not  calculated. 

The  second  change  implemented  as  a  result  of  the  Reference  7  data  base  has  to  do  with 
the  empirical  table  for  center  of  pressure  shifts.  Some  shift  changes  were  implemented,  which 
mainly  affect  results  in  the  transonic  region  for  lower  AOAs.  The  Reference  7  data  base  had 
Mach  0.9  data  available,  which  allowed  the  results  of  Reference  1  to  be  improved  upon 
somewhat.  These  modified  results  are  shown  in  Table  1.  They  result  in  some  slight 
improvement  in  the  average  center  of  pressure  error  for  the  Reference  1  data  base  from  about 
0.25  caliber  to  0.2  caliber.  The  0.2  caliber  error  is  an  average  error  of  about  1.6  percent  of  the 
body  length. 


The  third  body-alone  change  has  to  do  with  the  way  the  linear  and  nonlinear  terms  of 
Equation  (3)  are  treated  as  a  increases  above  30  deg.  The  AP98  methodology  assumes 


^n,  —  (Pn„  )ot,(X  —  30 


Cnl=(CklLo(i- 


cc-3  0^ 
60 


;30<a<90 


(6) 
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FIGURE  5.  COMPARISON  OF  MODIFIED  BODY-ALONE  AERODYNAMICS  METHOD 
TO  EXPERIMENT  FOR  FIGURE  3A  CASE 
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a  (deg.) 


FIGURE  6.  COMPARISON  OF  MODIFIED  BODY-ALONE  AERODYNAMICS  METHOD  TO 

EXPERIMENT  FOR  FIGURE  2A  CASE 
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TABLE  1  SHIFT  IN  BODY-ALONE  CENTER  OF  PRESSURE  AS  A  FUNCTION  OF  MACH  NUMBER 
AND  AOA  (AS  A  FRACTION  OF  BODY  LENGTH) 


a 

M 

0 

10 

20 

30 

40 

50 

60 

70 

80 

90 

0.00 

0.00 

0.025 

0.02 

0.000 

-0.025 

-0.040 

-0.040 

-0.030 

-0.010 

0.00 

0.20 

0.00 

0.025 

0.02 

0.005 

-0.025 

-0.040 

-0.045 

-0.030 

-0.010 

0.00 

0.40 

0.00 

0.03 

0.025 

0.005 

-0.025 

-0.040 

-0.050 

-0.030 

-0.015 

0.00 

0.60 

0.00 

0.03 

0.025 

0.00 

-0.035 

-0.055 

-0.070 

-0.050 

-0.030 

0.00 

0.80 

0.00 

0.030 

0.020 

-0.025 

-0.050 

-0.070 

-0.070 

-0.050 

-0.015 

0.00 

0.90 

0.00 

0.030 

0.020 

-0.02 

-0.050 

-0.070 

-0.070 

-0.040 

-0.015 

0.00 

1.00 

0.00 

0.02 

-0.01 

-0.02 

-0.040 

-0.040 

-0.040 

-0.030 

-0.005 

0.00 

1.15 

0.00 

0.02 

-0.01 

-0.02 

-0.020 

-0.025 

-0.030 

-0.025 

-0.005 

0.00 

1.30 

0.00 

0.02 

-0.01 

-0.01 

-0.010 

-0.010 

-0.010 

-0.005 

0.000 

0.00 

1.50 

0.00 

0.02 

-0.01 

0.000 

0.000 

0.000 

-0.010 

-0.005 

0.000 

0.00 

2.00 

0.00 

0.02 

0.02 

0.020 

0.015 

0.010 

0.005 

0.000 

0.000 

0.00 

2.50 

0.00 

0.02 

0.02 

0.02 

0.015 

0.010 

0.005 

0.000 

0.000 

0.00 

5.99 

0.00 

0.02 

0.02 

0.02 

0.015 

0.010 

0.005 

0.000 

0.000 

0.00 

^6.00 

0.00 

0.02 

0.02 

0.02 

0.015 

0.01 

0.005 

0.000 

0.000 

0.00 

As  seen  in  Figure  5,  which  shows  comparisons  up  to  40  deg  AOA,  Equation  (6)  yields  results 
that  are  slightly  low  compared  to  data  for  ct  >  30  deg.  In  reality,  the  linear  term  does  not  decay 
in  the  fashion  of  Equation  (6),  but  is  probably  more  parabolic  in  nature.  A  better  representation 
of  the  physics  is  therefore  assumed  to  be 


Cnl  -  (Pnb  )a’ a  -  30 
CNL=(cNLL30  30°<a<45= 

Cnl  =  (cNl  )a= 


(7) 


(,  oc -45 


1- 


45 


;  45  <  a  <  90° 


Equation  (7)  provides  some  additional  slight  improvements  in  body-alone  normal  force  and 
center  of  pressure  for  a  >  30  deg.  These  slight  improvements  are  shown  on  Figure  6  as  “AP02.” 
The  results  using  Equation  (6)  are  noted  as  “AP02  (Eq.  6).”  It  is  somewhat  hard  to  distinguish 
the  values  of  the  AP02  from  the  “AP02  (Eq.  6)”  in  Figure  6  because  of  the  scale.  As  an 
example,  for  M»  =  3.0,  a  =  40  using  Equation  (7)  (which  is  denoted  as  “AP02”),  Cn  =  9.38, 
whereas  when  Equation  (6)  is  used,  “AP02  (Eq.  6),”  CN  =  9.1.  These  values  compare  to  the 
experimental  data  value  of  9.4.  The  center  of  pressure  is  6.2  calibers  from  the  nose  tip  with 
Equation  (7)  and  6.3  calibers  with  Equation  (6).  The  experimental  data  is  6.2  calibers.  Similar 
improvements  of  using  Equation  (7)  versus  Equation  (6)  are  found  at  M»  =  4.5  of  Figure  6. 
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2.2  WING- ALONE  MODIFICATIONS 

The  wing-alone  methodology  of  Reference  1  assumed  the  wing-alone  normal  force  could 
be  predicted  from  a  fourth-order  equation  in  AOA.  That  is,  assuming  no  wing  camber, 

CNw  =a,aw  +a2<  +a3a3w  +  a4a^  (8A) 

a2  =  34.044  (CN)a=15.  -4.824  (CN)a=3S.  +0.426  (CN)o=60.  -6.412a,  (8B) 

a3  =  -88.240  (CN )a=15.  +  23.032  (CN  )a=35.  -  2.322  (CN  )a=60.  + 1 1 .464a,  (8C) 

a4  =53.219  (CN)a=I5.  -17.595  (CN)a=35.  +2.661  (CN)a=60.  -5.971a,  (8D) 


The  term  ai  of  Equation  (8)  is  the  value  of  wing-alone  lift  curve  slope  at  a  =  0  given  by  linear 
theory.  The  terms  (CN)a=15. ,  (CN)a=35„  and  (CN)o=60.  are  values  of  the  wing-alone  normal  force 

coefficients  at  a  =  15,  35,  and  60  deg,  respectively,  defined  by  the  data  bases  of  References  5,  6, 
and  14.  Above  aw  of  60  deg,  extrapolation  of  the  aerodynamics  at  aw  of  60  deg  is  used.  For 
more  details  of  the  method,  the  reader  is  referred  to  Reference  15. 

The  center  of  pressure  of  the  wing-alone  and  wing-body  normal  force  is  assumed  to  vary 
in  a  quadratic  fashion  between  its  linear  theory  value  near  cc  =  0  and  the  centroid  of  the  planform 
area  at  a  =  60  deg.  If  A  and  B  are  the  centers  of  pressure  of  the  linear  and  nonlinear  normal 
force  terms  (in  percent  of  mean  geometric  chord),  and  ccw  =  cl  +  6,  then  the  center  of  pressure  of 
the  wing-body  or  wing-alone  lift  is 


(XCT)™  =(Xcp)w  =  A  +  ^|aw|(B-A)+^oJw  (A-B)  (9) 

Equation  (9)  is  the  methodology  used  for  roll  position  of  0  deg.  For  roll  position  of  45  deg,  an 
equation  for  a  center  of  pressure  shift  was  derived  in  Reference  16  to  account  for  the  difference 
in  load  on  the  windward  and  leeward  planes.  This  shift  is  added  to  Equation  (9)  for  the  roll 
position  of  <J>=  45  deg  and  is 


(AXcp)^ 


r  + 

r  b  i 

cos(o)2  sin(2a) 

f  0.8a 

lCr+Cty 

U  3  J 

. 

l  65  J 

;  a  <  65  (10A) 


r  + 


b  } 

Cr 

c‘T 

U+ctJ 

u 

3  J. 

cosO2  sin(2a) ;  a  >  65° 


(10B) 


Equations  (10A)  and  (10B)  contain  a  correction  to  the  original  center  of  pressure  shift 
derived  in  Reference  16.  This  change  is  the  square  of  the  cos  (O)  term  in  Equation  (10),  whereas 
in  Reference  16,  the  cos  (<$)  term  was  to  the  first  power.  The  reason  for  the  square  is  the  fact 
that  the  cos  ($)  term  does  two  things.  First,  it  rotates  the  normal  force  to  a  plane  normal  to  the 


13 


NSWCDD/TR-99/1 16 


body  axis  as  opposed  to  being  normal  to  the  wing.  Second,  the  cos  ($)  term  rotates  the  radius 
vector  to  the  lateral  center  of  pressure  of  the  wing  from  the  <E>  roll  position  to  the  horizontal 
plane.  Reference  16  omitted  this  last  rotation,  causing  a  slightly  more  forward  center  of  pressure 
shift  at  roll  than  was  warranted.  As  already  mentioned,  one  of  the  keys  to  the  Reference  1 
method  was  the  development  of  the  wing-alone  normal  force  coefficient  tables  for  values  of  aw 
of  15, 35,  and  60  deg. 

The  NASA/MDAC7  wing-alone  data  base  had,  in  principle,  a  couple  of  advantages  over 
the  data  bases  used  to  develop  the  wing-alone  tables  at  a  =  15,  35,  and  60  deg  used  in  the  wing- 
alone  prediction  methodology  of  the  AP98.1  First  of  all,  the  Reference  7  data  base  measured 
wing-alone  data  for  a  =  0  to  90  deg  and  from  M.  =  0.6  to  4.0.  The  data  bases  comprising  the 
tables  in  Reference  1  consisted  of  several  different  sets  of  data  (see  References  5,  6,  and  14)  to 
cover  the  Mach  number  range  of  interest.  In  some  cases,  data  from  References  5,  6,  and  14  was 
available  only  to  60  deg  AOA,  and  in  some  data  bases,  the  data  tended  to  give  a  stall  effect  at 
higher  AOA  and  so  was  not  useable.  On  the  other  hand,  data  from  Reference  7  was  more  limited 
in  wing  planforms  considered  than  in  some  of  the  other  data  bases  (References  5, 6,  and  14). 

As  a  result  of  the  new  data  base  from  Reference  7,  it  was  decided  to  compare  the 
Reference  7  data  base  to  the  AP98  tables  as  well  as  the  Stallings  data,5  which  the  author  still 
believes  is  the  best  wing-alone  data  base  available.  Comparisons  were  made  as  a  function  of 
AOA,  aspect  ratio,  Mach  number,  and  taper  ratio.  Figures  7  and  8  compare  the  results  of  the 
Stallings  data  base5  and  the  recent  NASA/MDAC7  data  base  at  Mach  numbers  of  1.6  and  4.0, 
respectively,  for  fins  7  and  8  of  Reference  7.  Fin  7  is  of  aspect  ratio  1.0  with  taper  ratio  0,  and 
has  a  semispan  of  1.5  in.,  whereas  fin  8  is  of  aspect  ratio  2,  taper  ratio  0,  and  semispan  3.0  m. 
Also  shown  on  the  figures  are  the  results  from  the  AP98  method  and  revisions  to  the  wing-alone 
tables  to  be  incorporated  in  the  AP02.  Several  points  are  worthy  of  note.  First  of  all,  at  both 
M  =  1.6  and  M  =  4.0,  the  Reference  7  and  5  data  are  in  excellent  agreement  for  fin  8  up  to  AOA 
of  40  to  45  deg.  Above  a  =  45  deg,  the  Reference  7  data  stalls.  Also,  the  Reference  7  data  is 
consistently  about  10  percent  lower  than  the  Reference  5  data  for  fin  7  at  M  =  1.6  and  4.0.  It  is 
theorized  that  since  the  Reference  7  data  was  taken  with  a  splitter  plate  and  Reference  5  with  a 
sting,  the  differences  in  the  data  are  due  to  the  measurement.  It  is  suspected  that  for  the  lower 
semispan,  boundary  layer  buildup  ahead  of  the  fin  on  the  splitter  plate  is  the  source  of  the 
10  percent  lower  value  of  CNw  of  Reference  7  data  compared  to  Reference  5.  In  other  words, 

for  small  span  wings,  the  lower  dynamic  pressure  due  to  the  boundary  layer  near  the  root  chord 
has  more  of  an  effect  than  for  the  larger  span  wings.  This  effect  is  magnified  for  small  taper 
ratios  since  the  wing  cross-sectional  area  is  the  largest  at  the  root  chord.  It  is  not  known  why  the 
flow  stalls  above  about  45  deg  for  the  splitter  plate  results.  However,  this  was  the  case  for  most 
of  the  Reference  7  results.  As  a  result  of  these  two  phenomena,  it  was  decided  to  use 
considerable  judgement  before  using  any  of  the  Reference  7  results  for  the  1.5  in.  semispan  or 
for  any  span  above  a  =  45  deg.  The  final  point  to  be  made  in  viewing  Figures  7  and  8  is  that  the 
revised  values  of  CNw ,  which  will  be  incorporated  into  the  AP02,  are  closer  to  the  Stallings  data 

than  the  AP98.  The  AP98  had  intentionally  increased  the  values  of  CNw  somewhat  to  account 
for  the  fact  that  the  Stallings  data  was  taken  on  fairly  thick  wings  in  order  to  accommodate  many 
pressure  taps.  It  was  theorized  that  these  thick  wings  would  lower  CNw  unrealistically.  The 
revised  data  decreases  this  thickness  penalty  and  is  therefore  much  closer  to  the  Stallings  data. 
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FIGURE  7.  COMPARISON  OF  NASA/MDAC7  WING-ALONE  DATA  BASE  TO  THAT 

OF  STALLINGS5  (M  =  1.6) 
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FIGURE  8.  COMPARISON  OF  NASA/MDAC7  WING-ALONE  DATA  BASE  TO 
THAT  OF  STALLINGS5  (M  =  4.0) 
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After  analyzing  the  data  bases  of  References  5,  7,  and  14  through  plots  of  CNw  versus 
Mach  number,  CNw  versus  AOA,  and  CNw  versus  aspect  ratio,  general  ground  rules  were 

reached  on  how  to  revise  the  AP98  wing-alone  tables  based  on  the  new  Reference  7  data  base. 
These  ground  rules  are  listed  in  Table  2.  As  a  result  of  the  logic  in  Table  2,  only  slight  changes 
were  made  in  the  AP98  wing-alone  tables. 


TABLE  2.  GROUND  RULES  FOR  USING  NASA/MDAC7  WING-ALONE  DATA  TO  REVISE  AP981  TABLES 

1 .  Use  data  above  a  =  45  deg  with  caution,  because  of  inconsistency  with  Stallings  data  at  M  2  1.6  (caused  by  stall 
effect  in  Reference  7  data). 

2.  Ignore  data  on  two  smallest  fins;  boundary  layer  effects  apparently  decreased  CNw  compared  to  Stallings  data. 

3  If  Stallings  and  NASA/MDAC  data  are  both  available,  put  more  weight  on  Stallings  data  base. 

4.  At  M  =  1 .6,  Stallings  data  was  extrapolated  from  a  =  50  deg  to  a  =  60  deg.  Hence,  we  modified  data  slightly 
lower  in  most  cases.  However,  for  M  >  1.6,  data  was  available  to  a  =  60  deg. 

5.  Revise  X  =  1 .0  in  analogy  to  X  =  0  and  0.5  for  MDAC  data  and  previous  tables  since  X  =  1 .0  data  was  not 

available  from  NASA/MDAC  wind  tunnel  test. _ _ _ _ _ _ 


The  data  of  Reference  7  was  for  aspect  ratios  of  0.67, 1.33, 2.0,  3.0, 4.0,  and  6.0,  whereas 
that  of  Reference  5  was  for  aspect  ratios  of  0.5, 1.0, 2.0,  and  4.0. 

The  final  revised  set  of  wing-alone  tables  is  given  as  Tables  3  through  5.  These  tables 
will  be  a  part  of  the  AP02  and  will  replace  those  currently  used  in  the  AP98.  These  tables  are 
shown  plotted  in  Figures  9  through  11  for  taper  ratios  of  0,  0.5,  and  1.0  as  a  function  of  aspect 
ratio  and  Mach  number.  As  already  mentioned,  values  shown  here  are  close,  but  slightly 
different  than  those  used  for  the  AP98. 


2.3  REFINEMENTS  FOR  WING-BODY  AND  BODY-WING  INTERFERENCE  FACTOR 

NONLINEARITIES 

This  section  of  the  report  will  consider  refinements  in  the  empirical  factors  used  to  model 
the  nonlinearities  in  the  wing-body  and  body-wing  interference  factors  due  to  AOA.  No  changes 
will  be  made  in  the  nonlinear  empirical  constants  associated  with  the  interference  factors  due  to 
control  deflection,  since  the  Reference  7  data  base  did  not  consider  control  deflection^  as  a 
parameter.  Also,  the  focus  here  will  be  on  the  roll  orientation  of  <E>  =  0  deg  (fins  in  plus  “+”  roll 
orientation).  However,  when  changes  are  made  in  the  empirical  constants  for  $  =  0  deg,  the 
constants  for  $  =  45  deg  will  be  considered  for  change  in  a  complementary  way  to  the  $  =  0  deg 
results. 
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TABLE  3.  VALUES  OF  (cNw  )a=15<> 


ASPECT 

TAPER 

MACH  NUMBER 

RATIO 

RATIO 

0 

0.6 

0.8 

1.0 

1.2 

1.6 

2.0 

3.0 

4.5 

^6.0 

0.0 

.18 

.18 

.18 

.225 

.24 

.24 

.21 

.17 

.14 

.11 

*0.1 

0.5 

.19 

.19 

.19 

.225 

.24 

.24 

.21 

.17 

.14 

.11 

1.0 

.19 

.19 

.19 

.225 

.24 

.24 

.21 

.17 

.14 

.11 

0.0 

.28 

.29 

.30 

.32 

.32 

.24 

.18 

.16 

0.5 

0.5 

.39 

.41 

.415 

.43 

.43 

.45 

.38 

.30 

.22 

.19 

1.0 

.34 

.34 

.36 

.42 

.42 

.43 

.37 

.30 

.22 

.19 

0.0 

.43 

.44 

.46 

.50 

.54 

.46 

.42 

.32 

.22 

.18 

1.0 

0.5 

.47 

.50 

.55 

.65 

.66 

.58 

.45 

.34 

.24 

.21 

1.0 

.46 

.58 

.60 

.54 

.45 

.35 

.26 

.22 

0.0 

.55 

.59 

.65 

.72 

.70 

.62 

.50 

.34 

.27 

.23 

2.0 

0.5 

.56 

.59 

.66 

.76 

.80 

.68 

.54 

.40 

.30 

.27 

1.0 

.65 

.66 

.71 

.75 

.80 

.67 

.54 

.40 

.29 

.27 

0.0 

.65 

.66 

.71 

.79 

.83 

.70 

.59 

.39 

.31 

.26 

*4.0 

0.5 

.69 

.71 

.75 

.88 

.91 

.75 

.69 

.45 

.32 

.29 

1.0 

.69 

.71 

.75 

.88 

1 

1 

.67 

.45 

.31 

.29 
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TABLE  4.  VALUES  OF  (cNw  )a=JS, 


TABLE  5.  VALUES  OF  (CNw  )o=6Q0 
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FIGURE  11.  REVISED  WING-ALONE  DATA  BASE  ( X  =  1.0)  (Continued) 


To  better  understand  the  interference  lift  components,  it  is  instructive  to  examine  the  total 
normal  force  of  a  configuration  as  defined  by  Pitts  et  al.12  This  is  given  by 


CN  =<-'N„  +[(KW(B)  +KB(W))a+(kW(B)  +kB(W))5W](CNa)* 

+  [(Kt(B)  +KB(T))a  +  (kT(B)  +kB(T))5T](cNa)  +  CNt(v)  +<^i 
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The  first  term  in  Equation  (1 1)  is  the  normal  force  of  the  body  alone,  including  the  linear 
and  nonlinear  components;  the  second  term  is  the  contribution  of  the  wing  (or  canard),  including 
interference  effects  and  control  deflection;  the  third  term  is  the  contribution  of  the  tail,  including 
interference  effects  and  control  deflection;  and  the  last  terms  are  the  negative  downwash  effect 
on  the  tail  or  body  due  to  wing-shed  or  body-shed  vortices.  The  K’s  represent  the  interference  of 
the  configuration  with  respect  to  AO  A,  and  the  k’s  represent  the  interference  with  respect  to 
control  deflection.  Each  of  these  interference  factors  is  estimated  by  slender  body  or  linear 
theory.  As  such,  they  are  independent  of  AOA. 

The  terms  that  will  be  considered  in  this  report  for  refinements  are  Kw(B)>  KB(w),  Kt<b) 
and  KB(T).  These  four  interference  factors  are  defined  in  the  general  form 

K  =  Ksbt  +  AK(M„ ,  a,  AR,  X)  (12) 

LX 

The  first  term  of  Equation  (12)  is  defined  by  linear  theory  or  slender  body  theory,  whereas  the 
second  term  is  defined  by  utilizing  several  large  wind  tunnel  data  bases  to  back  out  the 
nonlinearities  as  a  function  of  Mach  number,  angle  of  attack,  aspect  ratio  and  taper  ratio.  The 
general  nonlinear  trend  of  those  two  interference  terms  is  shown  in  Figure  12.  This  general  trend 
is  basically  the  same  for  both  the  $  =  0  and  =  45  deg  roll  orientations.  However,  the  five 
tables  of  data  that  define  the  empirical  constants  in  the  equations  of  Figure  12  are  different  for 
$  =  0  and  <&  =  45  deg  for  both  KW(b>  and  KB(W).  As  already  discussed,  for  cruciform  missiles, 
SBT  gives  no  roll  independence  for  low  AOA  values  of  Kw(b>  and  KB(w> 

As  seen  in  Figure  12,  KW(B)  in  general  can  deviate  slightly  from  SBT  or  LT  near  AOA  of 
0  deg.  It  then  decreases  until  it  reaches  a  minimum  value  and  then  approaches  a  value  of  1.0  at 
high  AOA.  On  the  other  hand,  KB<w)  can  either  increase  or  decrease  past  AOA  of  0  deg. 
Eventually,  it  also  decreases  until  it  reaches  some  minimum  value  at  high  AOA.  The  physics  of 
what  occurs  in  this  nonlinear  behavior  and  the  details  of  the  interference  factor  nonlinearities  are 
given  in  References  13  and  16.  For  ease  of  reference  for  the  reader,  a  brief  discussion  of  the 
physics  of  the  flow  that  underlies  Figure  12  will  be  given  here. 

In  examining  the  nonlinear  models  for  Kwcb)  and  KB(w>  of  Figure  12,  it  is  instructive  to 
try  to  correlate  the  mathematical  models  with  the  physics  of  the  flow.  The  wing-body 
interference  factor  is  somewhat  easier  to  understand  than  the  body-wing  interference.  The  wing- 
body  experimental  data  show  that  at  low  Mach  number,  SBT  slightly  underpredicts  the 
experimental  data.  AS  AOA  is  increased,  KW(b>  starts  decreasing  and  in  some  cases  decreases 
below  its  wing-alone  value.  As  AOA  increases,  Kw(b>  approaches  its  wing-alone  value.  As 
Mach  number  increases,  the  positive  interference  lift  on  the  wing,  caused  by  the  presence  of  the 
body,  is  lost  faster  and  faster  as  AOA  increases.  That  is,  the  wing-alone  solution  is  recovered 
much  faster  at  high  Mach  number  as  AOA  increases,  than  at  low  Mach  number.  This  is  believed 
to  be  the  result  of  the  Newtonian  Impact  mechanism  where,  at  high  Mach  number,  the 
momentum  of  the  air  particle  is  lost  almost  entirely  upon  direct  impact  on  a  surface,  as  opposed 
to  wrapping  around  the  surface  and  carrying  some  of  the  momentum  with  it,  as  at  low  Mach 
numbers. 
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The  KB(W)  model  contains  body  vortex  effects,  nose-  and  wing-to-wing  shock  effects,  as 
well  as  the  usual  added  dynamic  pressure  of  the  body  caused  by  the  presence  of  the  wing.  While 
some  of  the  trends  in  Figure  12  can  be  rationalized,  others  cannot  except  in  light  of  these 
combination  effects.  The  alternative  to  simultaneously  modeling  several  physical  phenomena  is 
to  try  to  estimate  the  effects  of  the  body  vortices  and  nose-  and  wing-shock  interactions  and 
subtract  them  from  the  experimental  data  for  the  configuration  tested  and  then  add  them  back  in 
analytically  for  another  configuration  of  interest.  This  process  not  only  complicates  the 
methodology,  but  adds  additional  inherent  errors  because  these  effects  cannot  be  easily  and 
accurately  estimated.  The  present  approach  neglects  some  of  the  scale  effects  caused  by  the 
position  of  the  wing  on  the  body.  However,  this  error  is  probably  smaller  than  the  results  of 
approximating  analytically  the  other  effects,  subtracting  them,  and  then  adding  them  again  for  a 
different  geometrical  configuration. 

In  general,  KB(W)  actually  increases  with  AOA  at  low  Mach  numbers  to  a  certain  point, 
where  it  starts  decreasing  analogous  to  Kw(b>-  However,  a  certain  amount  of  lift  or  force 
enhancement  is  gained  all  the  way  to  a  =  90  deg  for  low  Mach  numbers  as  shown  in  Figure  12. 
This  phenomenon  is  assumed  to  occur  all  the  way  to  M  =  6.0  based  on  extrapolated  data  from  the 
point  where  experimental  data  end,  which  is  AOA  of  25  to  40  deg  depending  on  Mach  number, 
to  a  =  90  deg  shown  in  Figure  12. 

Additional  higher  AOA  data  above  a  =  40  deg  is  needed  for  both  KW(B)  and  KB(W)  to 
modify  the  assumed  extrapolations  of  the  models  for  Kw(B)  and  KB(W)  at  high  AOA.  However, 
until  additional  data  are  available,  the  approximate  nonlinear  models  for  KW(b>  and  KB(W>  can  be 
used  to  estimate  aerodynamics  for  engineering  use.  This  statement  will  be  validated  for  a  limited 
set  of  flight  conditions  in  a  later  section. 

The  way  the  nonlinearities  are  treated  for  the  second  term  of  Equation  (12)  is  by  using 
five  tables  for  AKw(B)  and  five  tables  for  AKB(W)  •  Also,  these  tables  are  different  for  0  =  0  and 
45  deg  roll  orientation.  These  tables  define  the  parameters  shown  in  Figure  12.  The  definition 
of  these  10  parameters  is  as  follows: 

[AKw(B)]cc=o  =  difference  between  SBT  and  data  at  a  =  0 
ac  =  angle  of  attack  where  Kw<b)  starts  decreasing 


d  KW(B) 
da 

CCD 

[AKB(w)]a=o 


=  rate  of  decrease  of  Kw(b>  between  a  =  ac  and  a  =  aD 

=  angle  of  attack  where  KW(b>  reaches  an  initial  minimum 
=  angle  of  attack  where  Kw(B)  reaches  a  constant  value 
=  difference  between  SBT/LT  and  data  at  a  =  0 
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^B(W) 

da 


ai 


«2 


[Kb(W)]m1N 


=  rate  of  change  of  KB(w>  between  a=  0  and  a  =  ai 

=  angle  of  attack  where  dKB(w/da  changes  sign 
=  angle  of  attack  where  KB(W)  reaches  a  constant 

=  constant  value  of  KB(w)  above  a  =  a2  as  a  percent  of  linear  theory  or  slender 
body  theory 


The  mathematical  models  for  Kw<b>  and  KB(w>  are  once  again  defined  based  on  SBT/LT 
and  the  empirical  data  for  the  constants  defined  previously.  The  specific  equations  for  Kw(B)  are 


^W(B)  —  [Kw(B)^BT  +  IAKw(B) 


(13  A) 


-  [^W(B)  L  +  lA^W(B)  +  I  (|a|  aC  )| 


dK 


W(B) 


da 


;  ac  <  a  <  ar 


(13B) 


1- 


f  cxM  ~H  ^ 

vaM-aD  j 


(l  [K-W(B)  ]a=otD  )>  aD  <  a  —  aM 


(13C) 


K-W(B)  -  [l^WCB)^^  ♦  a  >  aM  (13D) 

The  specific  mathematical  model  for  KB(w>  is  given  by  Equations  (14A)  through  (14C). 


For  a  £  ai, 


Kb(W)  —  ^(W)^  JaK b(W)  iotso  +H 


dK 


B(W) 


da 


(14A) 


For  a2, 


KB(W)  “  [^B(W)l,_ai  +  _a  { [^B(W) ]o:=0tl  [KbCW)^} 

i  ^2  1  J 


(14B) 


For  a  >  a2, 


KB(W)  -  [KB(W)Lin 


(14C) 
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Tables  6  through  15  give  the  revised  set  of  values  for  the  10  empirical  constants  of 
Figure  12  for  the  $  =  0  deg  roll  orientation  and  Tables  16  through  26  give  values  of  these  same 
constants  for  $  =  45  deg  roll. 

The  revised  values  of  the  empirical  constants  in  Tables  6  through  26  were  derived 
primarily  based  on  comparing  the  AP98  (including  the  revisions  of  Sections  2.1  and  2.2  of  this 
report)  to  the  wing-body  data  base  of  Reference  7.  The  empirical  constants  were  then  adjusted 
on  a  case-by-case  basis  to  improve  the  overall  predictions  of  theory  to  data.  Some  tables  were 
hardly  changed  from  those  of  Reference  1.  Other  tables,  such  as  Tables  12  and  23,  were  changed 
more  significantly. 


TABLE  6.  DATA  FOR  [AKw(B)la=o  -AT  $  —  0  DEG 


MACH  NUMBER 

ASPECT 

RATIO 

TAPER 

RATIO 

*0.1 

0.6 

0.8 

1.2 

1.5 

2.0 

2.5 

3.0 

3.5 

4.5 

*5.0 

*0.25 

0,0.5,  1.0 

0 

.25 

.25 

.15 

0 

0 

0 

0 

0 

0 

0 

0.5 

0.5 

.05 

.05 

.05 

.05 

0 

0 

0 

0 

0 

0 

0 

1.0 

0.5 

.25 

.15 

.05 

0 

0 

0 

0 

0 

0 

0 

0 

*2.0 

0.5 

.20 

.1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0.5 

0 

.30 

.35 

.2 

.18 

0 

0 

0 

0 

0 

0 

0 

1.0 

0 

.35 

.29 

.16 

.06 

0 

0 

0 

0 

0 

0 

0 

*2.0 

0 

.27 

.29 

.10 

.10 

0 

0 

0 

0 

0 

0 

0 

0.5 

1.0 

.05 

.05 

.05 

.05 

0 

0 

0 

0 

0 

0 

0 

1.0 

1.0 

.25 

.15 

.05 

0 

0 

0 

0 

0 

0 

o 

0 

*2.0 

1.0 

.20 

.1 

0 

.10 

0 

0 

0 

0 

0 

0 

0 
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TABLE  7.  DATA  FOR  ac  (deg)  AT  $  =  0  DEG 


MACH  NUMBER 

ASPECT 

RATIO 

TAPER 

RATIO 

<;0.1 

0.6 

0.8 

1.2 

1.5 

2.0 

2.5 

3.0 

3.5 

n 

*5.0 

*0.25 

0,  0.5, 1.0 

30.0 

22.0 

22.0 

10.0 

0 

0 

0 

0 

0 

0 

0 

0.5 

0.5 

30.0 

17.3 

11.5 

10.0 

0 

0 

0 

0 

0 

0 

0 

1.0 

0.5 

30.0 

20.0 

15.0 

10.0 

0 

0 

0 

0 

0 

0 

0 

*2.0 

0.5 

20.0 

15.0 

10.0 

15.0 

0 

0 

0 

0 

0 

0 

0 

0.5 

0 

20.0 

12.0 

10.0 

10.0 

0 

0 

0 

0 

0 

0 

1.0 

0 

40.0 

20.0 

15.0 

10.0 

0 

0 

0 

0 

0 

0 

0 

10.0 

20.0 

15.0 

15.0 

0 

0 

0 

0 

0 

0 

1.0 

30.0 

17.3 

10.0 

10.0 

0 

0 

0 

0 

0 

0 

0 

1.0 

1.0 

30.0 

15.0 

12.5 

10.0 

0 

0 

0 

0 

0 

0 

0 

*2.0 

1.0 

10.0 

15.0 

15.0 

15.0 

0 

0 

0 

0 

0 

0 

TABLE  8.  DATA  FOR  [r  w(B)  AT  $  —  0  DEG 


MACH  NUMBER 

ASPECT 

RATIO 

TAPER 

RATIO 

<;0.1 

0.6 

0.8 

1.2 

1.5 

2.0 

2.5 

3.0 

H 

*5.0 

*0.25 

0, 0.5, 1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

0.5 

0.5 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

.95 

1.0 

.97 

1.0 

1.0 

1.0 

0.5 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

*2.0 

0.5 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

0.5 

0 

1.0 

1.0 

1.0 

1.05 

.90 

.90 

.90 

.90 

1.0 

1.0 

0 

1.0 

1.0 

1.0 

.95 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

*2.0 

0 

1.0 

1.0 

.95 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

0.5 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.05 

1.15 

1.13 

1.15 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

.95 

.95 

.95 

1.0 

1.0 

1.0 

1.0 

1.0 

*2.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

.93 

.90 

.95 

1.0 
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TABLE  9.  DATA  FOR  aD  (deg)  AT  $  =  0  DEG 


MACH  NUMBER 

ASPECT 

RATIO 

TAPER 

RATIO 

<;0.1 

0.6 

0.8 

1.2 

1.5 

2.0 

2.5 

3.0 

3.5 

4.5 

£5.0 

*0.25 

0,0.5, 1.0 

80.0 

40.0 

38.0 

35.0 

35.0 

35.0 

30.0 

25.0 

20.0 

15.0 

15.0 

0.5 

0.5 

70.0 

33.0 

31.4 

27.5 

30.0 

16.8 

17.8 

17.0 

15.0 

15.0 

14.0 

1.0 

0.5 

60.0 

32.5 

44.0 

22.0 

20.0 

22.5 

17.5 

18.0 

10.0 

17.0 

15.0 

;>2.0 

0.5 

45.0 

35.0 

44.0 

40.0 

25.0 

16.5 

17.0 

16.0 

10.0 

12.0 

15.0 

0.5 

0 

25.0 

15.0 

14.0 

15.0 

15.0 

12.0 

11.5 

1.0 

0 

65.0 

31.0 

18.0 

21.5 

16.0 

17.0 

11.0 

13.0 

13.0 

£2.0 

0 

50.0 

35.0 

35.0 

30.0 

25.0 

20.0 

17.7 

17.0 

12.0 

12.6 

11.5 

0.5 

1.0 

70.0 

33.0 

34.2 

26.0 

30.0 

133 

13.4 

11.8 

12.2 

11.5 

1.0 

1.0 

60.0 

33.0 

40.0 

21.0 

20.0 

22.0 

17.0 

16.0 

9.0 

14.0 

12.0 

£2.0 

1.0 

45.0 

35.0 

35.0 

40.0 

25.0 

18.0 

15.0 

15.5 

12.0 

12.6 

11.5 

TABLE  10.  DATA  FOR  aM  (deg)  AT  $  =  0  DEG 


MACH  NUMBER 

ASPECT 

RATIO 

TAPER 

RATIO 

sO.l 

0.6 

0.8 

1.2 

1.5 

2.0 

2.5 

3.0 

3.5 

4.5 

£5.0 

*0.25 

0,0.5, 1.0 

80.0 

45.0 

45.0 

40.0 

44.0 

38.0 

50.0 

46.0 

50.0 

50.0 

46.0 

0.5 

0.5 

80.0 

33.0 

31.4 

40.0 

50.0 

17.0 

40.0 

17.0 

40.0 

15.0 

14.0 

1.0 

0.5 

80.0 

33.0 

45.0 

45.0 

50.0 

50.0 

50.0 

36.0 

33.0 

17.0 

17.0 

£2.0 

0.5 

80.0 

43.0 

45.0 

45.0 

50.0 

50.0 

50.0 

36.0 

33.0 

17.0 

17.0 

0.5 

0 

80.0 

30.0 

30.0 

40.0 

50.0 

48.0 

50.0 

50.0 

50.0 

50.0 

50.0 

1.0 

0 

80.0 

31.0 

40.0 

50.0 

50.0 

50.0 

44.0 

40.0 

40.0 

£2.0 

0 

80.0 

43.0 

45.0 

45.0 

50.0 

50.0 

50.0 

50.0 

50.0 

50.0 

35.0 

0.5 

1.0 

80.0 

33.0 

34.2 

50.0 

31.0 

50.0 

50.0 

50.0 

50.0 

50.0 

50.0 

1.0 

1.0 

80.0 

33.0 

40.0 

50.0 

42.0 

50.0 

50.0 

50.0 

44.0 

40.0 

40.0 

£2.0 

1.0 

80.0 

43.0 

45.0 

45.0 

25.0 

18.0 

15.0 

36.0 

33.0 

37.0 

30.0 

32 
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TABLE  11.  DATA  FOR  [AKB(W)]«=o  AT  4>  =  0  DEG 


MACH  NUMBER 

ASPECT 

RATIO 

TAPER 

RATIO 

sO.l 

0.6 

0.8 

1.2 

1.5 

2.0 

2.5 

3.0 

3.5 

H 

*5.0 

*0.25 

0,  0.5, 1.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.5 

0.5 

0.0 

-.28 

-.15 

.16 

.15 

.05 

0.0 

0.0 

0.0 

1.0 

0.5 

0.0 

-.20 

-.20 

.15 

.10 

.15 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

-.20 

-.07 

.1 

.18 

.10 

0.0 

0.0 

0.0 

0.0 

0.0 

0.5 

0 

0.0 

-.33 

-.30 

.28 

.20 

.10 

.08 

0.0 

0.0 

0.0 

0.0 

1.0 

0 

0.0 

-.24 

-.25 

.05 

.2 

.05 

0.0 

0.0 

0.0 

0.0 

0.0 

*2.0 

0 

0.0 

-.20 

0.0 

.17 

0.0 

0.0 

0.0 

0.0 

0.0 

0.5 

1.0 

0.0 

-.28 

-.15 

.13 

.15 

.10 

0.0 

0.0 

0.0 

0.0 

0.0 

1.0 

1.0 

0.0 

-.20 

-.20 

.22 

.10 

.05 

0.0 

0.0 

0.0 

0.0 

0.0 

*2.0 

1.0 

0.0 

-.20 

-.07 

.17 

.20 

.10 

.15 

0.0 

0.0 

0.0 

0.0 

TABLE  12.  DATA  FOR  dKB(w/da  (per  deg)  AT  <E>  =  0  DEG 


MACH  NUMBER 

ASPECT 

RATIO 

TAPER 

RATIO 

*0.1 

0.6 

0.8 

1.2 

1.5 

2.0 

2.5 

3.0 

3.5 

H 

*5.0 

*0.25 

0.0 

0.0 

0.0 

0.0 

-.006 

-.008 

-.010 

-.020 

-.024 

0.5 

0.5 

.003 

.023 

.023 

-.009 

-.018 

-.020 

-.015 

-.014 

-.015 

-.016 

-.020 

1.0 

0.5 

.003 

.012 

.006 

-.0075 

-.014 

-.016 

-.013 

-.014 

-.015 

-.020 

-.020 

*2.0 

0.5 

.006 

0.0 

0.0  . 

0.0 

-.008 

-.012 

-.014 

-.015 

-.016 

-.020 

0.5 

0 

.003 

.035 

.028 

0.0 

0.0 

0.0 

-.004 

i 

o 

£ 

-.015 

-.016 

-.020 

1.0 

0 

.003 

.020 

.0225 

-.0075 

-.011 

-.012 

-.013  ! 

-.014 

-.015  1 

-.020 

-.020 

*2.0 

0 

.003 

.008 

.006 

0.0 

0.0 

-.008 

-.012 

-.014 

-.015 

-.016 

-.020 

0.5 

1.0 

.003 

.038 

.033 

-.003 

-.010 

-.020 

-.015 

-.016  1 

-.020 

1.0 

1.0 

.003 

.007 

.005 

-.0075 

-.014 

-.016 

-.015 

-.016 

-.016 

-.020 

-.020 

*2.0 

1.0 

.003 

.006 

0.0 

■am 

-.015 

-.016 

-.020 

33 
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TABLE  13.  DATA  FOR  a,  (deg)  AT$  =  0  DEG 


MACH  NUMBER 

ASPECT 

RATIO 

TAPER 

RATIO 

*0.1 

0.6 

0.8 

1.2 

1.5 

2.0 

2.5 

3.0 

3.5 

H 

*5.0 

sO.25 

0, 0.5, 1.0 

15.0 

21.1 

16.5 

45.0 

23.3 

20.5 

18.0 

15.0 

10.0 

0.5 

0.5 

30.0 

22.2 

43.0 

40.0 

25.0 

25.0 

25.0 

20.0 

20.0 

0.5 

30.0 

25.0 

20.0 

70.0 

20.0 

0.00 

10.0 

10.0 

10.0 

10.0 

10.0 

*2.0 

0.5 

30.0 

20.0 

20.0 

40.0 

30.0 

30.0 

30.0 

24.0 

20.4 

26.0 

26.0 

0.5 

0 

30.0 

15.0 

15.0 

25.0 

25.0 

20.0 

20.0 

10.0 

27.0 

20.0 

20.0 

1.0 

0 

30.0 

25.0 

20.0 

70.0 

20.0 

0.00 

10.0 

10.0 

10.0 

10.0 

10.0 

*2.0 

0 

30.0 

25.0 

20.0 

40.0 

30.0 

30.0 

30.0 

32.0 

30.0 

20.0 

20.0 

0.5 

1.0 

30.0 

40.0 

25.0 

26.5 

21.6 

20.0 

20.0 

1.0 

1.0 

30.0 

25.0 

20.0 

70.0 

20.0 

0.00 

10.0 

10.0 

10.0 

10.0 

10.0 

*2.0 

1.0 

30.0 

ESI 

32.0 

26.0 

20.0 

26.0 

26.0 

TABLE  14.  DATA  FOR  a2  (deg)  AT  $  =  0  DEG 


MACH  NUMBER 

ASPECT 

RATIO 

TAPER 

RATIO 

sO.l 

0.6 

1.5 

2.0 

2.5 

3.0 

3.5 

H 

*5.0 

s0.25 

0,  0.5,  1.0 

90.0 

75.0 

65.0 

63.4 

60.0 

55.0 

52.5 

40.0 

47.5 

45.0 

42.5 

0.5 

0.5 

90.0 

75.0 

65.0 

62.0 

43.0 

41.0 

42.5 

25.0 

42.0 

40.0 

40.0 

1.0 

0.5 

90.0 

75.0 

75.0 

80.0 

40.0 

50.0 

40.0 

30.0 

30.0 

30.0 

30.0 

*2.0 

0.5 

90.0 

75.0 

75.0 

80.0 

90.0 

90.0 

42.0 

40.0 

40.0 

40.0 

40.0 

0.5 

0 

90.0 

75.0 

75.0 

80.0 

49.0 

47.8 

42.5 

43.0 

26.5 

40.0 

40.0 

1.0 

0 

90.0 

75.0 

75.0 

80.0 

40.0 

50.0 

40.0 

40.0 

30.0 

30.0 

30.0 

*2.0 

0 

90.0 

75.0 

75.0 

80.0 

90.0 

90.0 

41.0 

40.0 

40.0 

43.0 

43.0 

0.5 

1.0 

90.0 

75.0 

53.2 

41.0 

42.5 

26.5 

43.5 

40.0 

40.0 

1.0 

1.0 

90.0 

75.0 

74.0 

72.0 

40.0 

50.0 

40.0 

40.0 

30.0 

30.0 

30.0 

*2.0 

1.0 

90.0 

75.0 

75.0 

80.0 

lEH 

45.0 

30.0 

40.0 

43.0 

43.0 

34 
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TABLE  15.  DATA  FOR  [KB(W)]min  AS  A  FRACTION 
OF  SLENDER  BODY  THEORY  AT  $  =  0  DEG 


NL 

[Kb(W)]mIN 

0 

0.5 

3.8 

0.5 

4.9 

0.25 

6.0 

0 

TABLE  16.  DATA  FOR  [KW(b)]«=o  AT  $  =  45  DEG 


MACH  NUMBER 

ASPECT 

RATIO 

TAPER 

RATIO 

£0.1 

0.6 

0.8 

1.2 

1.5 

2.0 

2.5 

3.0 

3.5 

H 

25.0 

^0.25 

0,  0.5, 1.0 

0.00 

0.00 

0.00 

0.00 

0.00 

0 

0 

0 

0 

0 

0 

0.5 

0.5 

0.00 

0.00 

0.00 

-0.13 

0.00 

0 

0 

0 

0 

0 

0 

1.0 

0.5 

0.10 

0.00 

0.00 

0.00 

-0.10 

0 

0 

0 

0 

0 

0 

22.0 

0.5 

0.10 

0.00 

0.00 

0.00 

0.00 

0 

0 

0 

0 

0 

0 

0.5 

0 

0.00 

0.00 

0.00 

0.00 

0.00 

0 

0 

0 

0 

0 

:>2.0 

0 

0.00 

0.00 

0 

0 

0 

0 

0 

0 

0.5 

1.0 

0.00 

0.00 

0.00 

0.00 

0.00 

0 

0 

0 

0 

0 

0 

22.0 

1.0 

0.05 

0.00 

0.00 

0.00 

0.00 

0 

0 

0 

0 

0 

0 

1.0 

0 

0.35 

0.15 

0.05 

0.00 

-0.10 

0 

0 

0 

0 

0 

0 

1.0 

1.0 

0.10 

0.00 

0.00 

0.00 

-0.10 

0 

0 

0 

0 

0 

0 

35 


NSWCDD/TR-99/1 16 


TABLE  17.  DATA  FOR  ac  AT  3>  =  45  DEG 


MACH  NUMBER 

ASPECT 

RATIO 

TAPER 

RATIO 

*0.1 

0.6 

0.8 

1.2 

1.5 

2.0 

2.5 

3.0 

3.5 

H 

^5.0 

<;0.25 

0,0.5, 1.0 

0.0 

22.0 

22.0 

0 

0 

0 

0 

0 

0 

0.5 

0.5 

15.0 

11.5 

11.0 

10.0 

0.0 

0 

0 

0 

0 

0 

0 

1.0 

0.5 

15.0 

13.3 

0.0 

6.5 

0.0 

0 

0 

0 

0 

0 

0 

£2.0 

0.5 

10.0 

10.0 

0.0 

6.5 

2.2 

0 

0 

0 

0 

0 

0 

0.5 

0 

30.0 

15.0 

11.5 

10.0 

0.0 

0 

0 

0 

0 

0 

0 

£2.0 

0 

10.0 

10.0 

0.0 

6.5  ^ 

0.0 

0 

0 

0 

0 

0 

0 

0.5 

1.0 

15.0 

11.0 

11.0 

10.0 

0.0 

0 

0 

0 

0 

0 

0 

a  2.0 

1.0 

10.0 

10.0 

0.0 

6.5 

1.5 

0 

0 

0 

0 

0 

0 

1.0 

0 

EH 

0.0 

6.5 

0.0 

0 

0 

0 

0 

0 

0 

1.0 

1.0 

15.0 

13.3 

0.0 

6.5 

0.0 

0 

0 

0 

0 

0 

0 

TABLE  18.  DATA  FOR  AT  3>  =  45  DEG 


MACH  NUMBER 

ASPECT 

RATIO 

TAPER 

RATIO 

<;0.1 

0.6 

0.8 

1.2 

1.5 

2.0 

2.5 

3.0 

3.5 

H 

£6.0 

*0.25 

0,0.5,  1.0 

1.0 

1.0 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.0 

0.5 

0.5 

1.0 

1.0 

1.00 

0.90 

0.90 

1.00 

0.95 

1.00 

0.97 

1.00 

1.0 

1.0 

0.5 

1.0 

1.0 

1.00 

0.95 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.0 

£2.0 

0.5 

1.0 

1.0 

0.95 

0.95 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.0 

0.5 

0 

1.0 

1.0 

1.00 

1.00 

0.90 

0.90 

0.90 

0.90 

1.0 

£2.0 

0 

1.0 

1.0 

0.95 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.0 

0.5 

1.0 

1.0 

1.0 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.0 

£2.0 

1.0 

1.0 

1.0 

1.00 

1.00 

1.00 

1.00 

1.00 

0.93 

0.90 

0.95 

1.0 

1.0 

0 

1.0 

1.0 

1.00 

0.95 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.0 

1.0 

1.0 

1.0 

1.0 

1.00 

0.95 

0.95 

0.95 

1.00 

1.00 

1.00 

1.00 

1.0 
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TABLE  19.  DATA  FOR  aD  AT  0  =  45  DEG 


MACH  NUMBER 

ASPECT 

RATIO 

TAPER 

RATIO 

<;0.1 

0.6 

0.8 

1.2 

1.5 

2.0 

2.5 

3.0 

3.5 

4.5 

*6.0 

s0.25 

0,  0.5, 1.0 

20.0 

40.0 

38.0 

35.0 

30.0 

25.0 

16.3 

131 

13.1 

10.0 

0.5 

0.5 

59.0 

33.0 

30.0 

25.6 

25.0 

15.0 

15.0 

10.0 

15.0 

15.0 

10.0 

1.0  i 

0.5 

49.0 

38.0 

32.0 

26.0 

24.0 

14.4 

10.0 

10.0 

10.0 

*2.0 

0.5 

39.0 

31.5 

30.0 

28.0 

25.0 

16.5 

15.0 

14.4 

10.0 

13.0 

10.0 

0.5 

0 

59.0 

35.5 

33.0 

39.5 

29.5 

15.0 

25.0 

15.0 

15.0 

10.0 

10.0 

^2.0 

0 

39.0 

31.5 

30.0 

17.0 

13.5 

11.4 

10.0 

10.0 

10.0 

0.5 

1.0 

59.0 

35.5 

33.0 

25.6 

29.5 

15.0 

15.0 

15.0 

12.0 

13.0 

10.0 

22.0 

1.0 

39.0 

31.5 

30.0 

28.0 

23.3 

14.0 

16.0 

15.0 

11.8 

12.0 

10.0 

1.0 

0 

59.0 

38.5 

27.1 

17.2 

21.0 

11.4 

10.0 

10.0 

10.0 

1.0 

1.0 

49.0 

38.5 

32.5 

26.0 

26.4 

16.0 

15.3 

15.0 

11.8 

10.0 

10.0 

TABLE  20.  DATA  FOR  aM  AT  $  =  45  DEG 


MACH  NUMBER 

ASPECT 

RATIO 

TAPER 

RATIO 

sO.l 

0.6 

0.8 

1.2 

1.5 

2.0 

2.5 

3.0 

3.5 

H 

*6.0 

*0.25 

0, 0.5, 1.0 

35.0 

45.0 

45.0 

40.0 

44.0 

28.0 

25.0 

29.0 

20.0 

0.5 

0.5 

65.0 

33.0 

30.0 

49.0 

52.0 

40.0 

40.0 

30.0 

25.0 

25.0 

.20.0 

1.0 

0.5 

55.0 

38.0 

47.0 

49.5 

66.0 

48.5 

45.0 

41.0 

40.0 

10.0 

20.0 

*2.0 

45.0 

31.5 

40.0 

56.0 

57.0 

45.0 

45.0 

41.0 

40.0 

28.0 

20.0 

0.5 

0 

65.0 

35.5 

33.0 

50.0 

46.0 

30.0 

30.0 

50.0 

20.0 

*2.0 

0 

45.0 

31.5 

40.0 

56.0 

55.0 

58.5 

49.8 

44.2 

41.5 

28.5 

20.0 

0.5 

1.0 

65.0 

35.5 

33.0 

49.0 

52.0 

40.0 

28.0 

24.0 

21.0 

13.0 

20.0 

*2.0 

1.0 

45.0 

31.5 

40.0 

44.0 

32.0 

28.0 

20.0 

1.0 

0 

70.0 

38.5 

49.0 

63.0 

41.5 

40.0 

20.0 

1.0 

1.0 

55.0 

EH 

49.5 

60.0 

47.5 

40.0 

33.0 

32.0 

20.0 

20.0 
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TABLE  21.  DATA  FOR  [kw(b)  L=cxm  AT  $  =  45  DEG 


TABLE  22.  DATA  FOR  [AKB(w>]«=o  AT  $  =  45  DEG 
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TABLE  23.  DATA  FOR  dKB(W/da  (PER  DEG)  AT  $  =  45  DEG 
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TABLE  24.  DATA  FOR  ax  (DEG)  AT  $  =  45  DEG 


MACH  NUMBER 

ASPECT 

RATIO 

TAPER 

RATIO 

<;0.1 

0.6 

1.5 

2.0 

2.5 

3.0 

3.5 

H 

26.0 

*0.25 

0,0.5,  1.0 

10.0 

15.0 

15.0 

15.0 

10.0 

10.0 

10.0 

10.0 

10.0 

10.0 

10.0 

0.5 

0.5 

10.0 

57.0 

20.0 

23.0 

23.0 

15.0 

20.0 

15.0 

10.0 

10.0 

10.0 

1.0 

0.5 

10.0 

30.0 

15.0 

17.5 

15.0 

15.0 

15.0 

2  2.0 

0.5 

10.0 

15.0 

15.0 

15.0 

10.0 

10.0 

10.0 

12.0 

10.0 

10.0 

10.0 

0.5 

0 

23.0 

19.0 

20.0 

22.5 

15.0 

10.0 

10.0 

10.0 

0 

15.0 

15.0 

10.0 

10.0 

10.0 

12.0 

10.0 

10.0 

10.0 

24.0 

25.0 

25.0 

16.0 

20.0 

15.0 

10.0 

10.0 

10.0 

22.0 

1.0 

10.0 

15.0 

15.0 

15.0 

10.0 

10.0 

10.0 

12.0 

10.0 

10.0 

10.0 

1.0 

0 

20.0 

25.0 

30.0 

30.0 

15.0 

17.5 

15.0 

15.0 

15.0 

1.0 

1.0 

10.0 

10.0 

20.0 

25.0 

30.0 

30.0 

15.0 

17.5 

15.0 

15.0 

15.0 

TABLE  25.  DATA  FOR  a2  (DEG)  AT  $  =  45  DEG 


MACH  NUMBER 

ASPECT 

RATIO 

TAPER 

RATIO 

*0.1 

0.6 

0.8 

1.2 

1.5 

2.0 

2.5 

3.0 

3.5 

H 

26.O 

<;0.25 

0,0.5,  1.0 

35.0 

55.0 

50.0 

50.0 

45.0 

40.0 

30.0 

27.5 

25.0 

0.5 

0.5 

75.0 

65.0 

55.0 

44.0 

44.0 

36.0 

30.0 

20.0 

1.0 

0.5 

75.0 

35.0 

30.0 

30.0 

60.0 

60.0 

62.0 

80.0 

80.0 

80.0 

80.0 

0.5 

75.0 

65.0 

30.0 

30.0 

30.0 

60.0 

42.0 

45.0 

45.0 

0.5 

0 

75.0 

60.0 

50.0 

52.0 

40.0 

35.0 

44.0 

50.0 

36.0 

30.0 

20.0 

0 

30.0 

30.0 

30.0 

60.0 

62.0 

80.0 

42.0 

45.0 

45.0 

0.5 

1.0 

55.0 

42.0 

40.0 

38.0 

44.0 

40.0 

36.0 

30.0 

20.0 

22.0 

1.0 

75.0 

65.0 

30.0 

30.0 

30.0 

60.0 

62.0 

80.0 

42.0 

45.0 

45.0 

1.0 

0 

30.0 

30.0 

60.0 

60.0 

62.0 

80.0 

80.0 

80.0 

80.0 

1.0 

1.0 

30.0 

30.0 

60.0 

60.0 

62.0 

80.0 

80.0 

80.0 

80.0 
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TABLE  26.  DATA  FOR  [Kb(wj]min  (FRACTION  OF  SBT/LT)  AT  45  DEG 


MACH  NUMBER 

ASPECT 

RATIO 

TAPER 

RATIO 

sO.l 

0.6 

3.5 

^0.25 

0, 0.5,  1.0 

0.25 

0 

0.5 

0 

0 

IS 

3.0  COMPARISON  OF  MODIFIED  THEORY  TO  NASA/MDAC  DATA  BASE 


The  first  thing  one  likes  to  do  after  making  changes  to  the  APC  is  to  compare  the  new 
predictions  to  the  data  base  upon  which  the  changes  were  based  in  a  comprehensive  fashion. 
The  data  base  of  Reference  7  consists  of  Mach  numbers  0.6, 0.9,  1.2, 1.6,  2.0,  2.3,  2.96,  and  3.95 
for  the  twelve  fin  planforms  of  Figure  3B  mounted  on  the  body  of  Figure  3A.  Angles  of  attack 
from  0  to  20  deg  were  considered  at  the  subsonic  and  transonic  Mach  numbers,  whereas  AOAs 
to  30  deg  were  considered  at  the  supersonic  conditions.  Comparisons  of  the  new  predictions  to 
the  Reference  7  data  base  and  the  AP98  are  shown  in  Figures  13  through  24  for  each  of  the 
twelve  body-tail  cases  shown  in  Figure  3.  In  examining  Figures  13  through  24,  it  is  seen  that  the 
improvements  made  to  the  AP98  (listed  in  the  figures  as  the  AP02)  improve  the  normal  force 
coefficient  prediction  accuracy  compared  to  experiment  and  the  AP98.  In  a  quantitative  sense, 
the  errors  of  the  AP02  compared  to  experiment  were  measured  at  a  =  10,  15,  20,  25,  and  30  deg 
where  data  was  available.  The  error  here  is  defined  by 


Error(%)= 


^  Nop  ^  N-THEORY 


Cn 


x  100 


(15) 


These  errors  were  then  averaged  by  individual  Mach  number  and  for  all  Mach  numbers  for  the 
12  fins  of  Figures  13  through  24.  These  results  are  shown  in  Table  27.  As  seen  in  the  table,  the 
average  errors  on  normal  force  coefficient  prediction  are  less  than  10  percent  for  any  Mach 
number  and  under  5  percent  for  the  entire  data  base.  While  not  shown,  the  average  total  error  for 
the  AP98  on  normal  force  is  closer  to  7  percent.  While  this  is  still  under  the  quoted  average  error 
of  ±  10  percent,  it  is  considerably  higher  than  that  given  by  the  improvements  which  will  be  part 
of  the  AP02. 
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AR  =  2.0,  A  =  0,  r/s  =  0.5,  A^A,  2  .049 

s  NASA/MDAC  DATA 

- AP02 

- AP98 


FIGURE  13.  COMPARISON  OF  NASA/MDAC  WING-BODY  NORMAL  FORCE 
WITH  AP98  PREDICTIONS  (FIN  NO.  1) 
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AR  =  2.0,  X  =  0,  r/s  =  0.5,  A^A,  =  049 

©  NASA/MDAC  DATA 

- AP02 

- AP98 
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FIGURE  14.  COMPARISON  OF  NASA/MDAC  WING-BODY  NORMAL  FORCE 
WITH  AP98  PREDICTIONS  (FIN  NO.  2) 
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AR  =  4.0,  A.  =  0,  r/s  =  0.33,  A^A,  ~  .098 

©  NASA/MDAC  DATA 

- AP02 

- AP98 


FIGURE  14.  COMPARISON  OF  NASA/MDAC  WING-BODY  NORMAL  FORCE 
WITH  AP98  PREDICTIONS  (FIN  NO.  2)  (Continued) 
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FIGURE  15.  COMPARISON  OF  NASA/MDAC  WING-BODY  NORMAL  FORCE 
WITH  AP98  PREDICTIONS  (FIN  NO.  3) 
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AR  =  6.0,  A  =  0,  r/s  =  .25,  Aj/Ap  =  .147 

©  NASA/MDAC  DATA 

- AP02 

- AP98 


FIGURE  15.  COMPARISON  OF  NASA/MDAC  WING-BODY  NORMAL  FORCE 
WITH  AP98  PREDICTIONS  (FIN  NO.  3)  (Continued) 


NSWCDD/TR-99/1 16 


AR  =  1.33,  X  =  0.5,  r/s  =  0.5,  A2fl\  =  .074 

©  NASA/MDAC  DATA 

-  AP02 

- AP98 
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AR  =  1 .33,  X  -  0.5,  r/s  =  0.5,  Aj/Ap  =  .074 

©  NASA/MDAC  DATA 

-  AP02 

- AP98 


FIGURE  16.  COMPARISON  OF  NASA/MDAC  WING-BODY  NORMAL  FORCE 
WITH  AP98  PREDICTIONS  (FIN  NO.  4)  (Continued) 
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FIGURE  17.  COMPARISON  OF  NASA/MDAC  WING-BODY  NORMAL  FORCE 
WITH  AP98  PREDICTIONS  (FIN  NO.  5) 
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AR  =  4.0,  X  =  0.5,  r/s  =  .25,  Aj/Ap  »  .221 

©  NASA/MDAC  DATA 

- AP02 

- AP98 
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FIGURE  19.  COMPARISON  OF  NASA/MDAC  WING-BODY  NORMAL  FORCE 
WITH  AP98  PREDICTIONS  (FIN  NO.  7) 
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AR  =  0.67,  X  =  .5,  r/s  =  .5,  Aj/A,  «  .147 

©  NASA/MDAC  DATA 

- AP02 

- AP98 


FIGURE  22.  COMPARISON  OF  NASA/MDAC  WING-BODY  NORMAL  FORCE 
WITH  AP98  PREDICTIONS  (FIN  NO.  10)  (Continued) 
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AR  =  1 .33,  X  =  .5,  r/s  =  .33,  Aj/V  *  -295 

©  NASA/MDAC  DATA 

- AP02 

- AP98 


FIGURE  23.  COMPARISON  OF  NASA/MDAC  WING-BODY  NORMAL  FORCE 
WITH  AP98  PREDICTIONS  (FIN  NO.  11) 
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AR  =  1.33,  X  =  .5,  r/s  =  .33,  Aj/V  =  -295 

©  NASA/MDAC  DATA 

- AP02 

- AP98 


FIGURE  23.  COMPARISON  OF  NASA/MDAC  WING-BODY  NORMAL  FORCE 
WITH  AP98  PREDICTIONS  (FIN  NO.  1 1)  (Continued) 
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AR  =  2.0,  X  =  .5,  r/s  =  .25,  A2fln/Apton  *  .44 

©  NASA/MDAC  DATA 

- AP02 

- AP98 


FIGURE  24.  COMPARISON  OF  NASA/MDAC  WING-BODY  NORMAL  FORCE 
WITH  AP98  PREDICTIONS  (FIN  NO.  12)  (Continued) 
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TABLE  27.  AVERAGE  NORMAL  FORCE  ERRORS  OF  AP02 
COMPARED  TO  NASA/MDAC7  DATA  BASE  ($  =  0) 


MACH  NO. 

NO.  POINTS 

AVERAGE  ERROR  (PERCENT) 

0.6 

42 

7.2 

0.9 

42 

8.7 

1.2 

42 

3.7 

1.6 

60 

2.7 

2.0 

60 

3.0 

2.3 

60  1 

3.1 

2.96 

60 

3.7 

3.95 

60 

4.9 

TOTALS 

426 

4.4 

To  check  and  see  if  the  AP02  improvements  have  a  positive  or  negative  impact  on 
predictions  for  the  aerodynamics  of  the  NASA  Tri-Service  data  base,  Tables  28  and  29  were 
prepared.  The  Tri-Service  data  base  consisted  of  Mach  numbers  0.6,  0.8,  1.2,  1.5,  2.0,  2.5,  3.0, 
3.5,  and  4.5  with  AOA  up  to  25  to  40  deg,  depending  on  Mach  number.  Fins  and  body  tested  are 
shown  in  Figure  4.  The  highest  aspect  ratio  fins  were  very  small  (AR  =  4),  so  the  data  associated 
with  those  fins  was  not  considered  in  the  Tables  27  and  28  averages.  Also,  the  aspect  ratio  2.0 
fin  data  was  only  considered  above  M  =  1.5  in  the  averaging  process.  The  overall  average  error 
for  442  data  points  in  the  <&  =  0  deg  roll  is  3.4  percent.  This  compares  to  a  value  between  4  and 
5  percent  for  the  AP98.  The  $  =  45  deg  roll  results  are  presented  in  Table  29.  Here,  362  data 
points  were  considered  at  the  same  AOA  and  Mach  numbers  as  for  the  =  0  deg  roll  position. 
The  average  normal  force  error  for  each  Mach  number  is  less  than  10  percent  and  the  overall 
average  for  the  entire  data  base  is  3.5  percent.  Reference  16  shows  that  the  AP98  average 
accuracy  for  the  $  =  45  deg  roll  is  6.2  percent  for  Cn  and  1.2  percent  of  the  body  length  for 
center  of  pressure. 

Table  30  then  combines  the  results  for  Tables  27  through  29  into  an  overall  average. 
This  overall  average  error  is  less  than  4  percent,  with  the  worst  case  averages  being  in  subsonic 
and  transonic  flow,  where  matching  the  optimum  critical  crossflow  Reynolds  number  is  quite 
difficult.  In  scanning  over  the  1230  data  points,  it  was  seen  that  some  worse-case  errors  can 
approach  30  percent  in  the  subsonic  region,  even  when  we  try  to  utilize  the  best  crossflow 
Reynolds  number  for  body-alone  results.  The  flowfield  changes  when  wings  are  added,  so  the 
best  critical  crossflow  Reynolds  number  for  the  body  alone  may  be  different  than  the  optimum 
value  for  the  wing-body.  Generally  speaking,  the  worst-case  errors  at  supersonic  speeds  are  at 
low  angle  of  attack  where  experimental  data  corrections  for  nonzero  AOA  were  not  made. 
Errors  as  high  as  15  percent  were  seen.  However,  errors  of  this  magnitude  for  a  single  data  point 
were  quite  rare.  It  is  seen  that  the  improvements  based  on  the  Reference  7  data  base  earned  oyer 
to  the  Reference  4  data  base  as  well.  Hence,  the  overall  accuracy  of  the  AP02  in  predicting 
lifting  characteristics  of  missile  configurations  should  be  slightly  improved  over  the  AP98. 
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TABLE  28.  AVERAGE  NORMAL  FORCE  ERRORS  OF  AP02 
COMPARED  TO  NASA/TRI-SERVICE4  DATA  BASE  ($  =  0) 


MACH  NO. 

NO.  POINTS 

AVERAGE  ERROR  (PERCENT) 

0.6 

25 

3.2 

0.8 

30 

4.8 

1.2 

33 

3.6 

1.5 

63 

2.2 

2.0 

59 

3.5 

2.5 

58 

2.6 

3.0 

58 

3.7 

3.5 

57 

3.9 

4.5 

59 

3.6 

TOTALS 

442 

3.4 

TABLE  29.  AVERAGE  NORMAL  FORCE  ERRORS  OF  AP02  COMPARED 
TO  NASA/TRI-SERVICE4  DATA  BASE  ($  =  45  DEG) 


MACH  NO. 

NO.  POINTS 

AVERAGE  ERROR  (PERCENT) 

0.6 

22 

4.8 

0.8 

23 

7.5 

1.2 

27 

3.8 

1.5 

49 

3.0 

2.0 

49 

3.5 

2.5 

48 

2.5 

3.0 

49 

3.2 

3.5 

46 

3.7 

4.5 

49 

2.7 

TOTALS 

362 

3.5 
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TABLE  30.  AVERAGE  NORMAL  FORCE  ERRORS  OF  AP02 
COMPARED  TO  COMBINED  DATA  BASES4  7 


MACH  NO. 

NO.  POINTS 

AVERAGE  ERROR  (PERCENT) 

0.6 

89 

5.5 

0.8-0.9 

95 

7.2 

1.2 

102 

3.7 

1.5- 1.6 

172 

2.6 

2.0 

168 

3.3 

2.3-2.5 

166 

2.8 

2.96-3.0 

167 

3.6 

3.5-3.95 

163 

4.2 

4.5 

108 

3.2 

TOTALS 

1230 

3.8 

No  average  error  on  center  of  pressure  was  made  because  of  time  constraints.  However, 
suffice  it  to  say  that  the  average  center  of  pressure  error  for  the  AP98  on  the  NASA  Tri-Service 
data  base  was  less  than  2  percent  of  the  body  length.13  Improvements  made  in  normal  force 
should  only  improve  these  already  excellent  predictions.  Likewise,  no  improvements  in  axial 
force  were  sought,  as  we  were  quite  happy  with  the  power-off  predictions  of  axial  force  from  the 
AP98.  Improvements  in  power-on  axial  force  will  be  addressed  in  a  future  report. 


4.0  COMPARISON  OF  AP02  TO  CONFIGURATIONS  OUTSIDE  THE 
REFERENCES  4  AND  7  DATA  BASES 


While  the  average  accuracy  comparisons  of  CN  to  experiment  of  Tables  26  through  28  are 
impressive  for  a  semiempirical  code,  the  true  measure  of  success  is  based  on  the  ability  to 
accurately  predict  aerodynamics  on  a  wide  variety  of  configurations  outside  the  data  bases  upon 
which  the  empirical  nonlinearities  were  derived.  Several  cases  will  be  considered  over  a  variety 
of  flight  conditions  to  make  the  determination  of  whether  the  improvements  added  to  the  AP02 
are  generically  applicable  to  other  missile  configurations  and  if  they  improve  the  accuracy  of 
aerodynamic  estimation  over  the  AP98.  For  this  part  of  the  validation  effort,  roll  positions  of 
both  0  =  0  and  45  deg  will  be  considered.  Also,  comparisons  of  the  AP02  to  the  AP98  as  well 
as  to  wind  tunnel  data  will  be  given. 

The  first  case  considered  is  taken  from  Reference  17  and  is  a  model  of  an  older  version  of 
the  SEASPARROW  missile.  A  fairly  extensive  data  base  exists  for  this  configuration.  The 
configuration  is  shown  in  Figure  25,  where  the  wings  are  used  for  control. 
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FIGURE  25.  WING-BODY-TAIL  CONFIGURATION  USED  IN  VALIDATION  PROCESS 

(ALL  DIMENSIONS  IN  INCHES) 


This  configuration  has  a  length  of  about  18  calibers  with  a  tangent  ogive  nose 
2.25  calibers  in  length.  It  has  wings  and  tails  of  fairly  high  aspect  ratios  of  2.8  and  2.6 
respectively.  Data  was  taken  at  Mach  numbers  of  1.5  to  4.63  for  AOAs  to  40  deg  and  control 
deflections  of  0  and  10  deg  (at  M  of  1.5  and  2.0)  and  0  to  20  deg  (at  M  of  2.35  to  4.63).  The  data 
was  taken  at  a  Reynolds  number  of  2.5  x  106/ft  and  boundary  layer  trips  were  also  used.  The 
model  had  a  hollow  camber,  and  camber  axial  force  measurements  were  given  separately  in 
Reference  20.  These  results  were  added  to  the  forebody  axial  force  measurements  to  compare 
with  the  AP98  and  AP02. 

Figure  26  shows  the  comparisons  of  the  AP98  and  AP02  to  the  data  of  Reference  17  for 
0  =  0  deg  and  O  =  45  deg.  Figure  26A  and  26B  give  Ca,  Cn,  and  Cm  for  M„  =  1.5  at  6\v  =  0  and 
6W  =  10  deg  at  O  =  0  deg.  In  general,  both  the  AP98  and  AP02  give  good  comparisons  to  data. 
Figure  26C  and  26D  give  similar  results  for  M„  =  2.87,  and  Figure  26E  and  26F,  for  =  4.6. 
Overall,  for  this  configuration,  at  O  =  0  deg  roll,  the  AP02  and  AP98  are  about  equal  in  overall 
accuracy  comparisons.  The  worst  case  errors  are  for  center  of  pressure  at  higher  Mach  number 
and  AOA  where  the  bow  shock  intersects  the  wing  shocks.  This  nonlinear  phenomenon  is  not 
modeled  in  the  <3>  =  0  deg  roll  orientation  at  all.  For  the  0  =  45  deg  roll,  the  center  of  pressure 
shift,  Equation  (10),  partially  accounts  for  this  phenomenon,  but  not  entirely.  Center  of  pressure 
errors  approach  0.6  caliber  or  3  percent  of  the  body  length  at  M»  =  4.6  and  a  =  40  deg.  The 
other  point  is  that  the  normal  force  predicted  by  the  AP02  for  combined  oc  and  6\v  is  better  than 
the  AP98  at  M«  =  4.6,  whereas  the  opposite  is  true  for  the  axial  force  at  high  AOA.  The  reason 
for  this  phenomenon  has  to  do  with  the  fact  that  the  wing-alone  normal  force  tables  were 
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decreased  slightly  at  a  =  60  deg  for  the  AP02.  This  had  the  effect  of  lowering  the  configuration 
normal  force,  to  be  more  in  line  with  experimental  data,  but  moving  the  axial  force  at  combined 
a  and  Sw  away  from  the  data  (see  Figure  26F). 

Figures  26G  and  26H  present  M„  =  1.5  results  of  theory  compared  to  experiment  for 
O  =  45  deg  roll.  Figures  26G  and  26H  give  Ca,  Cn,  and  Cm  for  M  =  1.5  and  Sw  =  0  and  10  deg. 
Figures  261  and  26J  give  similar  results  for  M«  =  2.87  and  Sw  =  0  and  20  deg,  and  Figures  26K 
and  26L  give  results  for  =  4.6  and  Sw  =  0  and  20  deg.  Figure  26L  shows  the  same 
phenomena  as  Figure  26F:  that  is,  with  the  lowering  of  CNw  at  a  =  60  deg  in  Table  5  to  be  more 

in  line  with  the  Stallings  data,5  improvements  in  normal  force  coefficient  prediction  are  realized, 
but  at  the  expense  of  axial  force  coefficient  prediction.  Note  that  reasonably  good  agreement  is 
obtained  between  experimental  data  and  both  the  AP98  and  AP02  for  all  static  aerodynamics  at 
all  three  Mach  numbers  and  for  all  control  deflections.  Here,  the  worst-case  center  of  pressure 
error  is  less  than  3  percent  of  the  body  length. 

In  general,  for  the  configuration  of  Figure  25,  the  AP02  average  errors  show  only  slight 
improvements  over  the  AP98.  This  is  primarily  because  the  AP98  comparisons  to  data  were 
already  extremely  good  and  the  changes  to  the  AP98  methodology  based  on  the  Reference  7  data 
base  were  minor  for  this  configuration.  In  fact,  while  Cn  and  Cm  are  predicted  slightly  better 
with  the  AP02  than  with  the  AP98,  the  axial  force  at  combined  a  and  6W  is  slightly  worse  with 
the  AP02  than  with  the  AP98. 

The  second  configuration  is  taken  from  Reference  18  and  is  a  canard-body-tail  missile 
configuration.  It  is  22.2  calibers  in  length,  and  the  nose  is  hemispherical.  The  tail  surfaces  are 
fairly  large,  with  aspect  ratio  0.87,  and  fairly  thick,  with  truncated  trailing  edges.  The  canards 
are  aspect  ratio  1.73.  The  configuration  is  shown  in  Figure  27 A.  The  hangers  which  are  on  the 
wind  tunnel  model  were  not  modeled  by  the  APC.  Tests  were  conducted  for  M»  =  0.2  to  4.63, 
AOA  of  0  to  20  deg,  control  deflections  of  0  to  20  deg,  roll  of  0  to  45  deg,  Rn/ft  of  2  x  10  for  a 
model  with  boundary  layer  trips.  Base  pressure  values  as  a  function  of  M.  and  AOA  were  given 
in  Reference  18,  and  these  values  were  added  to  the  axial  force  information  so  total  axial  force 
values  could  be  shown. 

It  should  be  pointed  out  that  the  tail  thickness  in  Figure  27A  is  less  than  that  of 
Figure  32A  in  Reference  1.  Reference  1  incorrectly  used  the  value  of  0.236  in.  for  the  tail 
thickness,  versus  the  correct  value  of  0.109  in.  as  shown  in  Figure  27 A.  This  larger  value  of 
thickness  was  the  primary  source  of  the  overprediction  in  axial  force  coefficient  in  Reference  1 
using  the  AP98.  The  correct  value  of  tail  thickness  was  used  for  both  the  AP98  and  AP02 
computations  in  this  report. 

Figure  27B  gives  the  comparison  of  theory  and  experiment  for  3>  =  0  deg  roll  for  both  0 
and  20  deg  control  deflections.  Results  are  shown  in  terms  of  CA,  CN,  and  CM  versus  Mach 
number  for  a  =  20  deg.  Viewing  Figure  27B,  it  is  seen  that  the  AP98  and  AP02  both  give  good 
agreement  to  data.  In  comparing  the  AP02  to  the  AP98  and  experiment,  it  is  seen  that  the  AP02 
shows  some  improvement  in  prediction  of  normal  force  and  pitching  moment  coefficients 
compared  to  the  AP98  for  the  following  conditions:  (1)  Mach  numbers  less  than  0.9  and  (2) 
Mach  numbers  greater  than  2.1  for  normal  force  coefficient.  For  the  intermediate  Mach 
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FIGURE  26.  COMPARISON  OF  EXPERIMENT  AND  THEORY  FOR  CA)  CN  AND  CM 
FOR  FIGURE  25  WING  CONTROL  CASE  (Continued) 
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FIGURE  27A.  CANARD-BODY-TAIL  CONFIGURATION  WITH  HEMISPHERICAL  NOSE18 
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numbers,  prediction  accuracy  of  the  two  versions  of  the  APC  is  comparable.  Axial  force 
prediction  accuracy  for  this  configuration  of  the  two  codes  is  also  comparable,  since  only  minor 
changes  were  made  to  the  AP02  with  respect  to  axial  force  estimation. 

The  $  =  45  deg  roll  comparisons  of  Ca,  Cn,  and  Cm  for  a  =  20  deg  and  6c  —  0  and  20  deg 
are  shown  in  Figure  27C.  In  general,  the  AP02  gives  better  normal  force  coefficient  predictions 
compared  to  data  than  does  the  AP98.  Pitching  moment  coefficients  predicted  by  the  AP02  are 
also  slightly  better  than  those  predicted  by  the  AP98,  although  the  improvement  is  not  as  great  as 
for  the  normal  force  coefficient.  Again  little  difference  in  axial  force  coefficient  is  seen  between 
the  AP02  and  AP98. 

To  summarize  the  second  validation  case  considered,  it  is  seen  that  the  improvement  in 
normal  force  prediction  accuracy  of  the  AP02  based  on  the  more  recent  data  base  of  Reference  7 
carried  over  to  the  Figure  27 A  configuration.  For  the  56  data  points  of  Figures  27B  and  27 C  (14 
Mach  numbers,  2  roll  orientations,  and  2  control  deflections),  the  average  normal  force  error  was 
reduced  from  7.9  percent  using  the  AP98  to  4.2  percent  using  the  AP02.  This  is  a  reduction  in 
the  normal  force  prediction  error  of  over  40  percent.  Some  slight  improvement  in  pitching 
moment,  center  of  pressure,  and  axial  force  was  also  observed  for  the  AP02  compared  to  the 
AP98.  However,  these  improvements  were  not  nearly  as  large  as  for  normal  force  coefficient. 

The  third  configuration  was  tested  by  Jorgensen.19  The  configuration  is  shown  in 
Figure  28A  and  consists  of  a  wing-body  and  wing-body-tail.  Both  the  wings  and  tails  are  fairly 
large  in  surface  area  and  aspect  ratio.  Figure  28B  gives  the  normal  force  coefficient  comparison 
between  the  AP98,  the  AP02,  and  experiment  for  the  wing-body  case  at  Mach  numbers  of  0.6, 
0.9,  1.5,  and  2.0  and  AOA  to  60  deg.  The  AP02  provides  only  slight  overall  average  accuracy 
improvement  over  the  AP98.  Both  predictions  fall  well  below  the  average  accuracy  goal  of 
±  10  percent.  Figure  28C  gives  both  the  normal  force  and  center  of  pressure  comparisons  for  the 
wing-body-tail  case  of  Figure  28A.  Again,  the  AP02  shows  only  slight  improvement  over  the 
AP98  compared  to  experiment. 

The  fourth  case  considered  in  the  evaluation  of  the  improved  empirical  constants 
developed  for  the  nonlinear  aerodynamic  terms  of  the  normal  force  Equation  (1)  is  taken  from 
Reference  20  and  is  shown  in  Figure  29 A. 

The  wind  tunnel  model  was  about  22  calibers  in  length  with  a  sharp  nose  of  2.25  calibers. 
The  canards  had  an  aspect  and  taper  ratio  of  2.0  and  0.3  respectively.  Various  tail  fin  spans  were 
considered.  This  model  was  tested  at  Mach  numbers  1.6  to  3.5  at  AOA  to  about  18  to  20  deg.  It 
a  boundary  layer  trip  present  and  was  tested  at  a  RN/ft  of  2.0  x  106.  Reference  20  gave 
separate  values  of  base  axial  force  coefficient,  which  were  added  to  the  axial  force  values  given 
in  the  reference  to  compare  to  the  AP98  and  AP02  computations.  To  compare  the  experimental 
data  to  theory,  Mach  numbers  of  2.5  and  3.5  are  selected  at  roll  angle  45  deg.  Also,  values  of  the 
tail-to-canard  semispan  of  0.47  and  1.25  are  considered.  Figures  29B  and  29C  present  the 
comparison  of  theory  to  experiment  for  bt/bc  =  0.47  and  bt/bc  =  1.25  at  Mach  numbers  of  2.5  and 
3.5,  respectively,  for  Ca,  Cn,  and  Cm- 
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(AR)t  =  3.64  kT  =  29  d  =  2.6  in. 
(AR)W  =  2.81  Xt=:28 


FIGURE  28A.  WING-BODY  AND  WING-BODY-TAIL  CONFIGURATIONS 
USED  FOR  COMPARING  AP98  TO  EXPERIMENT  AND  AP02 
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Jorgensen  WB  Case  ($  =  0°) 


FIGURE  28B  NORMAL  FORCE  COEFFICIENT  COMPARISONS  FOR  WING-BODY 
CONFIGURATION  OF  FIGURE  28A 


82 


NSWCDD/TR-99/1 16 


FIGURE  29 A.  CANARD-BODY-TAIL  CONFIGURATION  WITH  VARYING  TAIL  SPAN 
(FROM  REFERENCE  20)  (ALL  DIMENSIONS  IN  INCHES) 


84 


NSWCDD/TR-99/1 16 


b/b,  =  0.47 


b«/bc  =  1.25 


1 — ^  j  9  ; 

|  j  |  j 

|  !  |  1  I 

i  |  |  !  : 

—  AP02 

'  : 

—  AP98 

1 . : . I . I . . . : . 

°  Exp.20 

2.5  5  7.5  10  12.5  15  17.5  20 

a 


7.5  10  12.5  15  17.5  20 
a 


2.5  5  7.5  10  12.5  15  17.5  20 


7.5  10  12.5  15  17.5  20 


10 

0 

-10 

-20 

-30 

-40 

-50 

-60 

-70 

,  -80 
-90 
-100 
-110 
-120 
-130 
-140 
-150 


5  7.5  10  12.5  15  17.5  20 

a 


0  2.5  5  7.5  10  12.5  15  17.5  20 

a 


FIGURE  29C.  COMPARISON  OF  THEORY  AND  EXPERIMENT  FOR 
CONFIGURATIONS  OF  FIGURE  29  A  ($  =  45  DEG,  M„  =  3.5) 
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In  examining  Figures  29B  and  29C,  it  is  seen  that  the  AP02  and  AP98  both  give  excellent 
agreement  with  experiment  for  the  bt/bc  =  0.47  case.  However,  for  the  bt/bc  —  1.25  case,  the 
AP02  shows  significant  improvement  over  the  AP98  in  both  Cn  and  Cm  at  both  M  =  2.5  and 
M  =  3.5.  Average  normal  force  coefficient  and  center  of  pressure  errors  were  reduced  by  a 
factor  of  two  or  more  for  this  case  with  the  AP02  compared  to  the  AP98  and  experiment. 

The  next  set  of  wind  tunnel  data  considered  for  comparison  purposes  is  taken  from 
Reference  21.  Body-alone,  body-tail,  and  wing-body-tail  configurations  were  all  a  part  of  this 
test  series.  Figure  30A  shows  one  of  the  configurations  tested  and  considered  here  for  validation 
of  the  AP02  results.  The  model  is  13.5  calibers  in  length  with  a  1.5  caliber  tangent  ogive  nose. 
The  wing  surfaces  are.  fairly  large,  with  thickness  of  t/cr  =  0.0178  and  wedge  angles  of  15  deg  on 
the  leading  and  trailing  edges.  The  tail  surfaces  have  thickness  of  t/cy  =  0.05  and  wedge  angles 
of  20  deg.  The  tests  were  conducted  at  Mach  numbers  0.7  to  3.08  with  Reynolds  number 
varying  from  about  2  x  106  to  4.6  x  106  per  foot.  The  smooth  model  without  boundary  layer  trip 
option  was  used  for  the  AP02  and  AP98  calculations.  AOA  to  25  deg  were  considered  in  the 
wind  tunnel  test.  For  comparison  purposes,  normal  force  and  pitching  moments  are  compared  to 
data  at  M„  =  1.42  and  3.08  for  the  <f>  =  0  deg  roll  orientation.  Figure  30B  presents  these  results. 
As  seen  in  the  figure,  both  the  AP02  and  AP98  give  quite  acceptable  comparisons  to  data,  with 
the  AP02  giving  slightly  better  comparisons  for  normal  force  coefficients.  The  pitching  moment 
prediction  of  the  AP02  and  AP98  are  about  equal  for  this  configuration  at  the  conditions 
considered,  with  the  AP02  being  slightly  better  at  M.  =  3.08  and  the  AP98  slightly  better  at 
M«o  =  1.42.  Both  versions  of  the  APC  give  aerodynamics  well  within  the  accuracy  goals. 
Reference  21  also  gives  axial  force  information  where  the  base  pressure  has  been  subtracted  out. 
Unfortunately,  only  a  side  camber  tap  was  used,  so  the  AOA  information  was  not  believed  to  be 
accurate.  Hence,  no  axial  force  comparisons  with  AOA  are  shown. 

Figure  31  shows  a  sixth  case  considered.  This  case  has  the  same  body  (12.33-caliber 
tangent  ogive-cylinder  with  a  3-caliber  nose)  as  that  tested  at  Langley.4  However,  dorsals  of 
aspect  ratio  0.1  and  tail  surfaces  of  aspect  ratio  2.0  have  been  added.  Mach  numbers  considered 
are  4.5  and  10.0.  This  case  was  originally  defined13  to  allow  computations  to  be  performed  with 
a  full  Euler  solver22  at  high  Mach  number,  since  wind-tunnel  data  above  Mach  4.6  appeared  to  be 
lacking.  Figure  3  IB  shows  the  comparison  of  the  ZEUS22  computations  for  normal  force  and 
center  of  pressure  with  the  AP98  and  AP02  results  for  the  body-tail  configuration,  and 
Figure  3 1C  shows  the  same  comparisons  for  the  body-dorsal-tail  case  of  Figure  31  A.  ZEUS 
computations  are  here  used  as  the  truth  model,  although  the  full  Euler  solution  also  has  some 
small  errors  in  comparison  to  wind  tunnel  data,  based  on  past  experience.  ZEUS  computations 
are  shown  only  for  cc  =  1, 10,  20,  and  30  deg.  It  is  seen  that  both  the  AP98  and  AP02  give  good 
comparisons  to  the  ZEUS  computations,  and  the  overall  average  comparisons  to  the  ZEUS 
computations  for  the  AP98  and  AP02  are  about  the  same. 

The  cases  considered  to  this  point  in  the  results  and  discussion  are  cases  that  lie  within 
the  parameter  space  of  the  Reference  7  data  base  except  for  the  Figure  31  configuration.  Since 
Tables  6  through  26  are  used  at  all  Mach  numbers  and  the  Reference  7  data  base  only  goes  down 
to  Moo  =  0.6,  it  was  decided  to  adjust  the  empirical  constants  in  the  tables  based  on  several 
configurations  outside  the  Reference  7  data  base  for  Mach  numbers  near  zero.  The  first  low 
Mach  number  case  is  shown  in  Figure  32 A  and  is  taken  from  Reference  23.  This  configuration 
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D  =  30  mm 


FIGURE  30A.  WING-BODY-TAIL  CONFIGURATION  CONSIDERED  FOR 
VALIDATION  WITH  AP02  AND  AP98  (REFERENCE  21) 
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FIGURE  32A.  MMPT  CONFIGURATION  TESTED  AT  M-  =  0.2  (FROM  SMITH,  SALAZAR  et  al.23) 


is  17.5  calibers  in  length  and  has  a  1.88-caliber  tangent-ogive  nose.  It  has  aspect  ratio  1.38 
wings  located  near  the  point  where  moments  are  taken,  and  aspect  ratio  1.56  tail  surfaces  located 
flush  with  the  base.  Axial  force,  normal  force,  and  pitching  moment  coefficients  predicted  by 
the  AP98  and  the  modified  AP98,  or  AP02,  are  compared  to  the  Reference  23  experimental  data 
in  Figure  32B.  Results  are  given  for  the  3>  =  0  deg  roll  plane  only.  Reference  23  also  stated  that 
axial  force  measurements  could  have  significant  errors  due  to  use  of  a  sting  designed  for  large 
normal  forces  at  high  AOA.  In  viewing  Figure  32B,  it  is  seen  that  both  the  AP02  and  AP98  give 
good  predictions  of  normal  force  and  pitching  moment  coefficients  compared  to  experiment. 
Axial  force  coefficient  predictions  for  the  AP02  and  AP98  are  identical  and  follow  the  trends  one 
would  expect,  although  discrepancies  with  experiment  exist  because  of  the  balance  used  for 
measurements.  In  a  quantitative  sense,  the  average  normal  force  coefficient  error  of  the  AP98 
was  7.4  percent  for  the  10  AOAs  considered.  The  AP02  reduced  this  average  error  to  5.0  percent 
or  about  a  one-third  reduction  in  average  error. 

The  next  low  Mach  number  case  is  shown  in  Figure  33,  and  the  test  data  was  given  in  a 
report  by  Howard  and  Dunn.24  This  configuration  has  dorsals  that  have  an  aspect  ratio  of  0.12 
and  tail  surfaces  that  have  an  aspect  ratio  of  4.0.  The  exact  configuration  illustrated  at  the  top  of 
Figure  33  is  not  within  the  allowable  constraints  for  fin  planform  required  by  the  APC. 
Therefore,  a  modified  version  of  the  fin  planforms  is  required,  one  that  meets  the  constraints  of 
the  APC.  This  configuration  is  shown  in  the  middle  of  Figure  33.  Note  that  the  parameters  that 
were  held  constant  for  the  fin  planforms  were  area,  aspect  ratio,  span,  taper  ratio,  leading-edge 
sweep  angle,  and  location  of  the  geometric  centroid  of  the  planform  area.  The  Howard  and 
Dunn24  work  gave  only  normal  force  as  a  function  of  AOA.  The  AP02  and  AP98  results  are  also 
shown  at  the  bottom  of  Figure  33.  Quite  acceptable  agreement  is  obtained  with  the  AP02 
compared  to  experiment,  even  at  high  AOA.  The  AP98  and  AP02  are  somewhat  lower  than  the 
data  suggest  at  high  a.  However,  part  of  this  underprediction  is  suspected  to  be  the  tendency  of  a 
base-mounted  sting  to  give  larger-than-true  normal  forces  at  subsonic  Mach  numbers.  In 
making  this  statement,  sting  interference  effects  were  assumed  to  be  unaccounted  for  in 
Reference  24.  Comparing  the  results  of  the  AP02  to  the  AP98  in  a  quantitative  sense,  the 
average  normal  force  error  of  the  AP98  for  34  data  points  is  10.7  percent,  whereas  the  average 
normal  force  error  of  the  AP02  is  6.0  percent.  This  6.0  percent  error  is  based  on  34  data  points  at 
both  the  $  =  0  and  <E>  =  45  deg  roll  orientations. 
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M  =  0.2 


FIGURE  32B.  COMPARISON  OF  STATIC  AERODYNAMIC  COEFFICIENTS  BETWEEN  EXPERIMENT 
AND  PREDICTIONS  FOR  FIGURE  32A  CONFIGURATION  ($>  =  0  DEG,  M.  =  0.2) 


94 


95 


NSWCDD/TR-99/1 16 


The  last  case  considered  for  the  low  Mach  number  validation  of  the  refined  nonlinear 
empirical  constants  of  Tables  6  through  26  is  shown  in  Figure  34A.  It  was  also  tested  in  the 
Naval  Postgraduate  School  wind  tunnel  in  Monterey,  California.  This  case  is  a  canard-controlled 
configuration  with  aspect  ratio  1.59  canard  and  0.9  tail  surfaces  (see  Figure  34 A).  The  body  is 
22.6  calibers  in  length.  Data  are  available  for  M»  =  0.2  up  to  45  deg  AOA  and  for  control 
deflections  of  the  canards  of  -20  deg,  0  deg,  and  +20  deg  for  both  $  =  0  and  45  deg  roll.  Results 
will  be  presented  here  for  only  the  0  deg  control  deflection  case  at  roll  of  both  0  and  45  deg. 
Figure  34B  presents  axial,  normal,  and  pitching  moment  coefficient  comparisons  of  the  AP02 
and  AP98  to  experimental  data27  for  $  =  0  deg  roll.  The  same  comment  applies  here  as  to  the 
previous  case  with  respect  to  wind  tunnel  accuracy  of  axial  force  measurements  using  a  sting 
balance  system  designed  for  measurement  of  high  AOA  normal  force  loads.  In  fact,  even  the 
normal  force  loads  at  low  AOA  are  suspect  since  both  the  0  and  5  deg  AOA  loads  are  negative. 
In  general,  the  AP02  has  slightly  improved  average  normal  force  predictions  compared  to 
experimental  data  between  AOA  10  deg  and  45  deg  and  the  AP98.  The  AP98  and  AP02  give 
about  equal  results  for  pitching  moment  coefficients.  The  $  =  45  deg  roll  results  are  shown  in 
Figure  34C.  Again,  the  AP02  gives  slightly  improved  results  for  normal  force  predictions 
compared  to  the  AP98,  with  pitching  moment  predictions  being  about  equal.  Quantitatively,  the 
AP02  and  AP98  give  average  normal  force  prediction  errors  of  about  6.8  percent  and  8.2  percent 
respectively  for  the  25  data  points  at  a  =  10  to  45  deg  and  roll  of  0  and  45  deg. 


5.0  SUMMARY  AND  CONCLUSIONS 


To  summarize,  the  nonlinear  empirical  constants  used  in  the  APC  to  predict  nonlinear 
normal  force  and  pitching  moments  on  missile  configurations  at  high  AOA  have  been  refined 
based  on  a  more  recent  missile-component,  wind-tunnel  data  base.7  In  comparing  the  new 
aerodynamic  predictions  of  the  revised  code  (AP02)  to  the  latest  released  version  of  the  APC 
(AP98)  the  following  conclusions  were  drawn: 

(1)  The  refined  nonlinear  empirical  coefficients  reduced  the  average  normal  force  error 
of  the  AP02  compared  to  the  AP98  for  the  NASA/MDAC7  data  base  by  over  a  third 
(7.0  percent  average  to  4.4  percent  average  error  based  on  426  data  points) 

(2)  In  comparing  the  new  AP02  to  the  AP98  for  the  older  NASA  Tri-Service  Data  Base,4 
it  was  seen  that  the  improvements  made  to  the  empirical  constants  also  gave 
improvements  in  accuracy  of  normal  force  coefficient  for  this  data  base  as  well. 
Average  normal  force  errors  were  reduced  from  4  to  5  percent  for  the  AP98  to 
3.4  percent  for  the  AP02.  This  also  represents  close  to  a  one-third  reduction  in 
average  normal  force  coefficient  errors. 
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Note:  All  Dimensions  in  Inches,  Full  Scale 
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FIGURE  34A.  CANARD-CONTROLLED  MISSILE  CONFIGURATION27 
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FIGURE  34C.  COMPARISON  OF  AXIAL,  NORMAL  AND  PITCHING  MOMENT  COEFFICIENTS 
BETWEEN  EXPERIMENT,  AP02,  AND  AP98  FOR  FIGURE  34A  CONFIGURATION  ($  =  45  DEG,  M.  =  0.2) 
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(3)  No  quantitative  assessment  was  made  of  center  of  pressure  (or  pitching  moment) 
improvements.  However,  in  viewing  the  results  qualitatively,  it  is  believed  a  slight 
overall  improvement  was  realized  by  the  improved  normal  force  loads.  In  addition, 
an  error  in  the  center  of  pressure  shift  at  roll  of  45  deg  was  corrected  in  the  AP02, 
also  adding  some  slight  improvement  in  center  of  pressure  predictions. 

(4)  In  comparing  the  AP02  to  the  AP98  on  nine  wing-body-tail  configurations  outside  of 
the  missile  component  data  bases  upon  which  the  nonlinear  empirical  constants  were 
derived,  it  was  found  that  in  general,  the  improvements  in  average  normal  force  error 
of  the  AP02  were  seen  here  as  well.  The  average  improvements  range  from  only  a 
slight  improvement  on  one  case  to  over  40  percent  reduction  in  error  for  the  best  case. 
Overall,  it  is  guessed  that  the  average  normal  force  error  was  reduced  by  about  20  to 
25  percent  from  the  AP98  to  the  AP02. 

(5)  While  the  overall  accuracy  improvement  in  normal  force  coefficient  is  based  on 
averages,  one  can  still  find  a  single  data  point  error  on  either  the  AP98  or  AP02 
where  the  error  is  as  high  as  35  percent.  These  worst-case  data  points  usually  occur 
at  subsonic  or  transonic  speeds  where  it  is  very  difficult  to  predict  the  correct  value  of 
critical  crossflow  Reynolds  number  and  Mach  number. 

(6)  No  assessment  of  axial  force  or  control  deflection  errors  were  made  since  no  changes 
were  implemented  into  the  AP02  compared  to  the  AP98. 

(7)  Based  on  the  overall  improvement  in  normal  force  using  the  refined  nonlinear 
constants  of  this  report,  these  improvements  will  be  a  part  of  the  next  version  of  the 
aeroprediction  code  that  will  be  transitioned  to  users  in  Fiscal  Year  2002. 
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AOA 

APC 

AP02,  AP98 

AR 

LT 

NASA/LRC 

NASA/MDAC 


NSWCDD 
SB,  SBT 
Aref 


Aw 

b 


CA 

cAb.cAf,cAw 


C 


Cd 


Cl 


7.0  SYMBOLS  AND  DEFINITIONS 

Angle  of  Attack 
Aeroprediction  Code 

2002  and  1998  versions  of  the  APC  respectively 
Aspect  Ratio  =  b2/Aw 
Linear  Theory 

National  Aeronautics  and  Space  Administration/Langley  Research  Center 

National  Aeronautics  and  Space  Administration/McDonnell  Douglas 
Corporation 

Naval  Surface  Warfare  Center,  Dahlgren  Division 
Slender  Body,  Slender-Body  Theory 

Reference  area  (maximum  cross-sectional  area  of  body,  if  a  body  is 
present,  or  planform  area  of  wing,  if  wing  alone)(ft ) 

Planform  area  of  wing  in  crossflow  plane  (ft ) 

Wing  span  (not  including  body)(ft) 

Axial  force  coefficient 

Base,  skin-friction,  and  wave  components,  respectively,  of  axial  force 
coefficient 

Crossflow  drag  coefficient 

Drag  coefficient 
Lift  coefficient 
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CM 


'Ml 


'Mn 


Cn 

CM 


NB(V) 


CN  ,  CN 

NB(W)  WB(T) 

Cnl 

CNnl 

fc'Na)w’(pNet)T 

c. 


^T(V) 


CNw 

r  c 

NW(B)  ’  NF(B) 

cl 


Ct 


cal 


dref 


Pitching  moment  coefficient  (based  on  reference  area  and  body  diameter, 
if  body  present,  or  mean  aerodynamic  chord,  if  wing  alone) 

Linear  component  of  pitching  moment  coefficient 
Nonlinear  component  of  pitching  moment  coefficient 

Normal  force  coefficient 

Normal  force  coefficient  of  body  alone 

Negative  afterbody  normal-force  coefficient  due  to  canard  or  wing-shed 
vortices 

Normal-force  coefficient  on  body  in  presence  of  wing  or  tail 

Linear  component  of  normal-force  coefficient 

Nonlinear  component  of  normal-force  coefficient 

Normal  force  coefficient  slope  of  wing  and  tail  respectively 

Negative  normal-force  coefficient  component  on  tail  due  to  wing  or 
canard-shed  vortex 

Normal  force  coefficient  of  wing  alone 

Normal-force  coefficient  of  wing  or  fin  in  presence  of  body 

Normal-force  coefficient  derivative 

Root  chord  (ft) 

Tip  chord  (ft) 

Caliber(s)  (one  body  diameter) 

Body  diameter  (ft)  at  base 
Reference  body  diameter  (ft) 
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dKw(Ei- ,  dKe(Wl  Rate  at  which  KW(B)  or  KB(W)  decreases 
da  da 


deg  Degree(s) 

KB(w>,  KB(T)  Ratio  of  additional  body  normal-force  coefficient  in  presence  of  wing,  or 

tail-to-wing  or  tail-alone  normal-force  coefficient  at  6  =  0  deg 

kB(W),  kB(D  Ratio  of  additional  body  normal-force  coefficient  due  to  presence  of  wing 

or  tail  at  a  control  deflection  to  that  of  wing  or  tail  alone  at  a  =  0  deg 


[KB(W)]min  Minimum  value  of  KB(W)  as  percent  of  slender-body  theory  value 

KW(b),  KT(b>  Ratio  of  normal-force  coefficient  of  wing  or  tail  in  presence  of  body  to 

that  of  wing  or  tail  alone  at  6  =  0  deg 

kW(B),  kf(B)  Ratio  of  wing  or  tail  normal-force  coefficient  in  presence  of  body  due  to  a 

control  deflection  to  that  of  wing  or  tail  alone  at  a  =  0  deg 


AK 

[AKw(B)]ct=0 

and 

[AKB(W)]o=o 

f,  «n 


MN 


M, 


Nonlinear  component  of  wing-body  or  body-wing  interference 

Amount  that  the  experimental  values  of  Kw<b)  and  Kb(w>  exceed  slender 
body  theory  at  a  =  0  deg 

Body  length  and  nose  length  respectively 
Mach  number  normal  to  body  =  M»  sin  a 

Normal  Mach  number  where  flow  transitions  from  subcritical  to 
supercritical  conditions 

Freestream  Mach  number 


N  Normal  force 

NF  Fin  normal  force 


r 


Local  body  radius  (ft) 


RBM  Root  bending  moment 

rn  Reynolds  number  where  flow  transitions  from  subcritical  to  supercritical 

conditions 
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rw,  rj 


s 

XCP 


(Xcp)l.  (Xcp)nl 

x,y,z 

a 


«c 


aD 


aw,  cct 


«1,  <*2 
6 


6w,  6t 


*1 


TlO 

$ 


A 


Radius  of  body  at  wing  or  tail  locations 

Wing  or  tail  semispan  plus  the  body  radius  in  wing-body  lift  methodology 

Center  of  pressure  (in  feet  or  calibers  from  some  reference  point  that  can 
be  specified)  in  x  direction 

Center  of  pressure  of  linear  and  nonlinear  terms  of  normal  force 
Axis  system  fixed  with  x  along  centerline  of  body 
Angle  of  attack  (deg) 

Angle  of  attack  where  wing-body  interference  factor  starts  decreasing 
(deg) 

Angle  of  attack  where  the  wing-body  interference  factor  reaches  a 
minimum  (deg) 

Angle  of  attack  where  Kw<b)  reaches  a  constant  value 

Local  angle  of  attack  of  wing  or  tail  (aw  +  6  or  ax  +  6,  respectively,  in 
degrees) 

Angles  of  attack  used  in  nonlinear  model  for  Kb(W) 

Control  deflection  (deg),  positive  leading  edge  up 

Deflection  of  wing  or  tail  surfaces  (deg),  positive  leading  edge  up 

Parameter  used  in  viscous  crossflow  theory  for  nonlinear  body  normal 
force  (in  this  context,  it  is  the  normal  force  of  a  circular  cylinder  of  given 
length-to-diameter  ratio  to  that  of  a  cylinder  of  infinite  length) 

Value  of  t)  at  Mn  =  0 

Roll  position  of  missile  fins  (<J>  =  0  deg  corresponds  to  fins  in  the  plus  (+) 
orientation;  $  =  45  deg  corresponds  to  fins  rolled  to  the  cross  (x) 
orientation) 

Taper  ratio  of  a  lifting  surface  =  ct/cr 
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Subscripts 
c,  t,  w 
CG 


Canard,  tail,  wing 
Center  of  gravity 
Freestream  conditions 
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